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Abstract 


With the increase in the complexity of power system, mostly in the generation, 
control and transmission levels, stable operation of the system under various conditions is 
of main concern. Power system is highly non-linear and unpredictable in nature. So 
stability studies of the system as a whole is very challenging. Proper modeling of dynamics 
of all the components with sound knowledge of all the constraints is necessary for a 
detailed study of power system. With reasonable assumptions and simplifications the study 
can be made easier. 

In many cases, instability and eventual loss of synchronism are initiated by some 
spurious disturbance in the system resulting in oscillatory behaviour that, if not damped, 
may eventually build up. This is very much a function of the operating condition of the 
power system, such as strength of the system, generator power output, and Automatic 
Voltage Regulator (AVR) settings. Oscillations, even if undamped at low frequencies, are 
undesirable because they limit power transfers on transmission lines, in some cases, may 
take the generator out of synchronism and induce stress in the mechanical shaft. Power 
System Stabilizer (PSS) aids in damping these small oscillations via modulation of the 
generator exciter. The art and science of PSS application has developed considerably over 
last few decades. 

In this thesis the small signal instability problem associated with Kiaga Nuclear 
Power Plant (NPP) is addressed. The causes for the oscillation are determined and PSS is 
tuned to damp out these oscillations. An intelligent Fuzzy Logic PSS is also designed, 
which can be used in place of the existing PSS for proper operation of the system. Other 
standard methods are also studied and developed for application purposes. 



Chapter 1 


Introduction 


1.1 General 

Nuclear power generation is identified as one of the major electric power source in 
future. These power plants are very sensitive to any disturbance (mechanical and electrical) 
due to various technical and safety reasons. If the exciter and Power System Stabilizer 
(PSS), any, are not tuned properly, there may be sustained oscillations, which is very 
dangerous. These oscillations are related to power system stability. 

Power system stability is the defined as that property of power system that enables it 
to remain in a state of operating equilibrium under normal operating conditions and to 
regain an acceptable state of equilibrium after being subjected to a disturbance. In general 
stability is a condition of equilibrium between opposing forces. Instability results when a 
disturbance leads to a sustained imbalance between the opposing forces. In the evaluation 
of stability the concern is the behaviour of the power system when subjected to a transient 
disturbance. The disturbances may be small or large. Small disturbances in the form of load 
changes take place continuously, and the system adjusts itself to the changing conditions. 
The system must be able to operate satisfactorily under these conditions and successfully 
supply the maximum load. It must also be capable of surviving numerous disturbances of a 
severe nature, such as short circuit on a transmission line, loss of a large generator or load, 
or loss of a tie between two subsystems. The understanding of stability problems is greatly 
facilitated by the classification of stability into various categories. 

The classification can be based upon physical nature/main system parameters such as 
rotor angle, voltage and frequency. Based on the type of disturbances and time span for 
which a particular disturbance lasts, it can be further subdivided into various categories as 
shown in Fig. 1.1. 
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1.2 Rotor Angle Stability 


Rotor angle stability is the ability of interconnected synchronous machines to remain 
in synchronism after being subjected to a disturbance. It depends on the ability to restore 
equilibrium between electromagnetic torque and mechanical torque of each synchronous 
machine. If the system is perturbed such as change in load or input, any kind of fault in the 
system [21], this equilibrium is disturbed, resulting in acceleration or deceleration of the 
rotors of the machines according to the laws of motion of a rotating body. If one generator 
temporarily runs faster than another, the angular position of its rotor relative to that of the 
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slower machine will advance. The resulting additional angular difference transfers part of 
the load from the slow machine to the fast machine, depending on the power angle 
relationship [20]. This tends to reduce the speed difference and hence the angular 
separation. Beyond a certain limit, an increase in angular separation is accompanied by 
decrease in power transfer; this increases the angular separation further and leads to 
instability. 

For any given situation, the stability of the system depends on whether or not the 
deviations in angular positions of the rotors results in sufficient restoring torques. When a 
synchronous machine losses synchronism with the rest of the system, its rotor runs at a 
higher or lower speed than that required to generate voltage at system frequency. Loss of 
synchronism can occur between one machine and the rest of the system or between groups 
of machines. In the later case synchronism may be maintained within each group after its 
separation from others. The change in electrical torque of a synchronous machine following 
a perturbation can be expressed into two components: 

AT C = K^AS + jKj^Aco (11) 

where, KsAS is the component of torque change in phase with rotor angle perturbation AS 
and is referred to as the synchronizing torque component. K$ is the synchronizing torque 
coefficient or stiffness of the system. TqAco is the component or torque in phase with the 
peed deviation Aco and is referred to as the damping torque component; T D is the damping 
torque coefficient. 

Stability depends on the existence of both components of torque for each of 
synchronous machines. Lack of sufficient synchronizing torque results in instability 
through an aperiodic drift in rotor angle. On the other hand, lack of sufficient damping 
torque results in oscillatory instability. For convenience in analysis the rotor angle stability 
phenomenon can be divided into two categories: 

1.2.1 Small Signal (or small-disturbance) Stability 

It is the ability of the power system to maintain synchronism under small 
disturbances. Such disturbances occur continuously on the system because of small 
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variations in loads and generation. The disturbances are considered sufficiently small for 
linearization of the system equations to be permissible for purpose of analysis. Instability 
that may result can be of two forms: (1) steady increase in rotor angle due to lack of 
sufficient synchronizing torque, or (2) rotor oscillations of increasing amplitude due to lack 
of sufficient damping torque. The nature of system response to small disturbances depends 
on a number of factors including the initial operating condition, the transmission system 
strength and the type of excitation control used. With high gain and continuously acting 
voltage regulators, the small signal stability problem is one of the ensuring sufficient 
damping of the system. The high gain regulators improve the steady state stability limit but 
reduce the damping torque of the system [19]. The stability of the following types of 
oscillation is of concern 

1 . Local Modes are associated with the swinging of units at a generating station 
with respect to the rest of the power system. 

2. Inter Area Modes are associated with the swinging of many machines in one 
part of the system against machines in another part. They are caused by two or 
more groups of closely coupled machines being connected by weak ties. 

3. Control modes are associated with generating units and other controls. Poorly 
tuned exciters, speed governors, HVDC converters and static var compensators 
are the usual causes instability of these modes. 

4. Torsional Modes are associated with the turbine-governor shaft system 
rotational components. Instability of torsional modes may be caused by 
interaction with excitation controls, speed governors, HVDC controls, and series 
capacitor compensated system. 

1.2.2 Transient Stability 

It is the ability of the power system to maintain synchronism when subjected to a 
severe transient disturbance. The resulting system response involves large excursions of 
generator rotor angles and is influenced by the non-linear power angle relationship. 
Stability depends on both the initial operating state of the system and the severity of the 
disturbance. Usually the system is altered so that the post-disturbance steady state operation 
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differs from that prior to the disturbance. The system is however designed and operated so 
as to be stable for a selected set of contingencies. 

1.3 Voltage Stability and Frequency Stability 

It is the ability of power system to maintain steady acceptable voltage at all buses in 
the system under normal operating conditions and after being subjected to a disturbance. A 
system enters into a state of voltage instability when a disturbance, increase in load 
demand, or change in system conditions causes a progressive and uncontrollable drop in 
voltage. The main factor causing the instability is inability of the power system to meet the 
demand for reactive power. For purposes of analysis the voltage stability can be classified 
into two categories, (1) Large disturbance voltage Stability, and (2) Small disturbance 
voltage stability. Large disturbance voltage stability is concerned with a system’s ability to 
control voltages following large disturbances such as system faults, loss of generation, or 
circuit contingencies. Small disturbance voltage stability is concerned with a system’s 
ability to control voltages following small perturbations such as incremental change in 
system load. 

Frequency stability is the ability to maintain steady frequency within a nominal 
range following a disturbance resulting in a significant imbalance between generation and 
load. Instability that may result occurs in the form of sustained frequency swings leading to 
tripping of generating units and/or loads. Depends on the ability to restore balance between 
generation and load of island systems with minimum loss of load and generation. 
Generally, frequency stability problems are associated with inadequacies in equipment 
responses, poor coordination of control and protection systems. At present voltage stability 
and small signal oscillations have appeared as main problem. Increasing steady state 
stability by high AVR gain is affecting badly the small signal stability. 

1.4 State of The Art 

With increase in demand and beginning of power markets, the major thrust of interest 
is in the area of reliable operation of power system. Advancement in VLSI technology and 
computational facilities has encouraged the utilities to change the old technology. Almost 
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all the present day generating units are equipped with fast acting voltage regulators. 
Demello and Concordia [1] have studied the phenomenon of stability as affected by 
thyristor type excitation system. The basics regarding application of the stabilizing signal to 
damp the small signal oscillations are given by E.V. Larsen and D.A. Swann [2], The signal 
requirements, the potential problems like torsional interaction and noise are addressed in 
the three part paper [2]. Practical approach to supplementary stabilizing signal from 
accelerating power is given in [3]. Some alternative method to take care of the mechanical 
power variations during disturbance are given in [3], In reference [4], power system 
stabilizer is applied for enhancement of over all system stability, i.e. not only small signal 
stability, but also the transient or large signal stability. 

All the methods discussed in above references address the problem of PSS tuning 
based on a Single Machine Infinite Bus (SMIB) analysis, where the generator suffering 
form oscillations is considered to be connected to a very large system having very high 
inertia constant. But in practical situations, the controllers of different machines interact 
with each other during disturbance conditions. So this dynamic interaction has to be 
considered while tuning the PSS for a particular machine. This leads to analyze the whole 
stability problem for a multimachine system, where the dynamics of different generators 
are considered simultaneously. Multimachine PSS tuning considering output feedback and 
decentralization constraints are studied in [5]-[6]. In reference [7]-[8], the LQR problem is 
solved for placing the under damped eigen values in a vertical strip. Use of tie line power 
How as input for a decentralized controller to stabilize a longitudinal power system is 
studied in [9]. In 1992 IEEE reviewed the earlier standard for exciter models used for 
stability studies [10]. This standard also gives detail about the PSS structures used in 
practice. 

With advance in control system, the Conventional PSS (CPSS) can be replaced by 
intelligent PSS, which can automatically change its parameters as per the situation. In 
reference [14] Malik etal have coordinated CPSS and Fuzzy Logic PSS (FLPSS). 
Application of Neural network to prepare the rule base, to tune the FLPSS parameters are 
given in [13]. In references [15]-[16] Genetic Algorithm (GA) has been applied to optimize 
the PSS input and output parameters for improvement in system performance. 
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Implementation of a FLPSS is given in reference [17]. In reference [18], coordination of 
CPSS and Flexible AC Transmission Systems (FACTS) controller parameters are studied. 
References [19]-[21] have described the aspect of power system stability in detail. 
Reference [19] provides ideas regarding modeling of different devices such as generator, 
exciter, governor system, induction motor, and transformer and transmission lines for 
stability studies. This also describes PSS tuning for small signal oscillation problems, 
transient stability analysis, voltage stability analysis and torsional interactions. Reference 
[20] deals with dynamic stability and control of power systems. The detail modeling of 
DAE equations describing power system dynamics for both small signal and transient 
stability analysis is given in [21]. This book also describes the energy function approach for 
solution of stability problems. In reference [22], fuzzy logic is applied to power system 
operations, such as Load frequency scheduling, PSS tuning, optimal load flow etc. The 
practical implementation of Fuzzy logic for all these applications is given in detail. 
Reference [26] describes application of fuzzy logic to other fields, such as pattern 
recognition, level control etc. 

1.5 Motivation 

Kaiga Nuclear Power Plant (NPP) under Nuclear Power Corporation of India Ltd. 
(NPCIL) has an installed capacity of 2X220 MWe. It is connected to two nearby hydro 
generators at Kadra and Kudasalli by a single line operating at 220kV, as given in Fig. 1.2. 

It is also connected to Guttur 400kV substation with a long line of length 165 km length. 
One Inter Connecting Transformer (ICT) 220/400 kV is installed at Kaiga and Kaiga- 
Guttur line operates at 400kV. Earlier this line was operating at 220kV. Both the generators 
of Kaiga NPP have potential source feed type static exciter [15]. For faster response and to 
provide higher ceiling voltage during faults, the Automatic Voltage Regulator (AVR) has 
PID characteristics. The AVR gain and time constant can be varied over wide range of 
values by selector switches on the field. The generator, exciter, current setting of AVR 
parameters and other data are given in Appendix-C. 


7 




Fig. 1 .2: Kaiga NPP and nearby generators 

The Kaiga generators have experienced severe power oscillations, whenever both 
generating units are operating and 220kV Kaiga- Kadra and Kaiga-Kudasalli lines are not 
in operation and power was evacuated through the DC Kaiga-Guttur 220 kV line. The 
range of observed frequency of oscillation was between 0.8-1. 2 Hz. The oscillations, 
therefore, appear to be in local plant mode. No oscillation was observed in nearby Hydro 
generators. This indicates that generators at Kaiga are oscillating against a large power 
system connected through weak transmission line. In view of this it is inferred that PSS on 
the generators at Kaiga need adjustment to avoid such oscillation. An attempt has been 
made in this thesis to find out the different reasons for oscillations in the Kaiga generators. 
Later the existing PSS, which is not functional at present, is tuned to aid in damping these 
oscillations. All the system data, model and diagram are taken from [27] provided by the 
manufacturer BHEL. 
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1.6 Thesis Organization 


This thesis has been organized into four parts. 

❖ Chapter- 1 gives a brief idea about system stability. A state of the art in the area 
of the research work related to stability is described. The problem with Kaiga 
NPP is defined which is the motivation behind the present work. 

❖ Chapter-2 discusses the models of different major components of a power 
system. The differential algebraic equations of the power system are linearized 
for small signal analysis. Different PSS tuning methods have been described and 
applied to WSCC 9-bus system and New England 39-bus system. 

❖ Chapter-3 discusses about the existing AVR and PSS structures at Kaiga NPP. 
Linearized model of AVR is used along with other DAE equations for eigen 
value analysis of Kaiga system. Phase and gain plot are found out by an SMIB 
analysis method. Then PSS is tuned and eigen values were calculated for 
compensated system. 

❖ Chapter-4 introduces a new area of PSS design by application of intelligent 
control system, such as Fuzzy Logic. The basics of FLC design is described in 
brief. The rule base for different sets of inputs is prescribed. Finally the FLPSS 
is tuned and system response is observed. 

1.7 Paper 

1 . “ Power System Stabilizer Tuning in Nuclear Power Station: A Case Study of India ”, 
under preparation for IEEE conference. 
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Chapter 2 


Power System Modeling for PSS Tuning 

2.1 Introduction 

The designing of a large interconnected system to ensure stable operation at minimum 
cost is a very complex problem. From a control theory point of view, a power system is a 
very high-order multivariable process, operating in a constantly changing environment. 
Because of high dimensionality and complexity of the system, it is essential to make 
simplifying assumptions and to analyze specific problems using the right degree of detail of 
system representation. This requires a good grasp of the characteristics of the overall 
system as well as of those of the individual devices. The power system is highly nonlinear 
and stochastic in nature whose dynamic performance is influenced by a wide array of 
devices with different response rates and characteristics. Virtually the dynamic behaviour 
of every major element of the power system has an effect on the system stability. 
Knowledge of these characteristics is essential for the understanding and study of power 
system stability. In this chapter, modeling of power system components and power system 
stabilizers tuning are presented. 

2.2 Modeling of Power System Components 

For the dynamic study of a system, suitable dynamic model of all the devices should 
be considered. In case of power system, the dynamics of the system is highly nonlinear 
because of devices like dynamic loads (mainly induction motor loads), and it is very 
difficult to predict the system behaviour under all the operating conditions due to varying 
complexity of the system. In most of the situations, for a change in system conditions, some 
devices take part to their fullest extent where as some devices have a very little contribution 
or no contribution at all. Considering this aspect of dynamics we can take these major 
equipments only with some assumptions to simplify the analysis. 
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2.2.1 Synchronous Machine 


In most of the studies, balanced, symmetrical, three-phase synchronous machine with 
a field winding and two damper windings on the rotor (2-Axis model or 2.1 model [20]) is 
considered. Fig. 2.1. shows two-axis (fourth order) model of a synchronous machine with 
stator transients neglected [21]. 



Fig. 2. 1 : Synchronous machine two-axis model dynamic circuit 


The differential equation (in pu) describing the dynamics of m generators of a system 
is given as: 



( 2 . 1 ) 


da, = Zk [E' v -X'*I a \ [g'd-JTWj D M-coJ 

dt M, M, qi M, di M, 


Wd, 

T doi 



dE 'di 

dt 


E'u, -X\‘)l q , 

T' q „i T'coi 


where, 


( 2 . 2 ) 

(2.3) 

(2.4) 
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<5 = Machine rotor angle with respect to a synchronously rotating frame in degrees 

co = Rotor speed in rad/sec 

M=Machine inertia in sec/electrical radian 

T,„ = Mechanical input torque 

E q — Voltage across d-axis due to field flux 

E d = Voltage across q-axis due to damper winding flux across d-axis 

I q = Quadrature axis component of stator current 

Id = Direct axis component component of stator current 

X q = Quadrature axis stator steady state reactance 

Xd - Direct axis stator steady state reactance 

X d~ Direct axis stator transient reactance 

X q =Quadrature axis stator transient reactance 

T do ~ Direct axis open circuit time constant 

T qo = Quadrature axis open circuit time constant 

2.2.2 Exciter Modeling 

The physical device that provides the field voltage (Ef<i) is called exciter. The main 
function of exciter is to provide sufficient DC power to the field of the synchronous 
generator for production of field flux linkage. Output of the exciter is nothing but the field 
voltage of the synchronous machine but proper scaling should be done as the power rating 
of the two is different [19]. If a synchronous machine is to be useful for a wide range of 
operating conditions, it should be capable of maintaining constant voltage and frequency. 
Mechanical input of a machine controls the frequency. For long-term stability studies the 
dynamics of governor system is to be considered. If it is specified as constant field voltage 
and mechanical torque input, then it implies that machine does not have voltage and speed 
(frequency) control. 

In steady state, adjustment of field voltage changes the machine terminal voltage. As 
the field circuit time constant is high (of the order of few seconds), fast control of the field 
current requires field forcing. Thus exciter should have a high ceiling voltage, which 
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enables it to operate instantly with voltage levels that are 4 to 5 times the normal value. 
Depending on the method of power supply to the exciter, the exciters can be categorized as: 

■ DC excitation system which uses DC generator with brushes and slip rings, driven 
by main shaft. 

* AC excitation system, which uses either an inverted AC generator (brushless) 
driven off the main shaft with rotating rectifier or an AC generator with a static 
rectifier to produce DC power. 

■ Static exciter, such as AC-DC converter in which power is supplied from the 
synchronous machine or auxiliary power. 

Efd normally is not controlled directly, but it is changed through the actuation of the 
exciter or pilot exciter. The main exciter may have a pilot exciter that provides the means 
for changing the output of the main exciter. In order to automatically control terminal 
voltage, a transducer signal must be compared to a reference voltage and amplified to 
produce the exciter input (Vr). This amplifier comprises the Automatic Voltage Regulator 
(AVR). The amplifier can be a pilot exciter (another DC generator) or a solid state control. 
If the voltage input to the AVR amplifier is simply the error voltage produced by the 
difference between a reference voltage and a conditioned potential transformer connected 
to the synchronous machine terminals, the closed loop control system may exhibit 
instabilities. In any standard excitation system this is accomplished through a stabilizing 
transformer whose input is connected to the output of the exciter and whose output is 
subtracted from the amplifier input. This feed back stabilization is a part of AVR loop. In 
case of rotating type of exciters (the first two of the above classification) the relationship 
between the exciter voltage and current is non-linear because of saturation of the exciter 
iron. This is taken into account by a proper saturation function which varies as the field 
voltage Efd represented as Se (Efd) as shown in Fig. 2.2. This function is mainly taken as 
exponential function [19]. 

In literature the most commonly used exciter is IEEE - 1 [21] exciter. It essentially 
represents a DC exciter, with some modifications can also be used for representing static 
exciter. The regulator output Vr is limited which also imposes limit on Efd. Usually the later 
are usually specified and the former can be found from the dynamic equations applying the 
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steady state conditions. The dynamic equations of IEEE type-1 exciter of i th generator are 
given as: 



Fig. 2.2: IEEE Type- 1 exciter model 


dE iL= | y u 

dt r H , 7), 


dV 

Z_Jji 

dt 


V,, K 




At 


T T 

1 Ai 1 


R, 


Ai 




Em + 


K 


Ai 


T 

1 Ai 


<y«rV,) 


dt 


T» (T,f 


where, 

Tn = Exciter time constant in Second 
Ke = Exciter gain 
Vr = Regulator output, 

Maximum value of regulator output in PU 
Vr "" 1 = Minimum value of regulator output in PU 
T a = Regulator amplifier time constant 
K a = Regulator gain 
Rj= Rate feedback in PU 
7> = Feedback circuit time constant, sec 
Ky = Feedback circuit gain 


(2.5) 

( 2 . 6 ) 

(2.7) 
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2.3 Complete Differential-Algebraic Equations 


The two-axis model equations are rewritten after neglecting the sub-transient reactances 
and saturation. The turbine and governor dynamics are also neglected resulting in constant 
Tm,. The limit constraints on Vr, (AVR output) are also deleted, since the objective is to 
simplify die system for modeling and simulation. A linear damping torque is assumed, that 
is, Tfwi ~ A (coi -cos). The resulting Differential Algebraic Equations (DAE Eqs) for an m- 
machine, n-bus system with the IEEE-type I exciters are given as follows: 


1 . Differential Equations for i th machine is 


dS, 

dt ~ C ° i C ° s 


dco, 

T u , [E‘ 

qi - X dil di ) j 



D,(co, -co,) 

dt 

M, 

M, 


M, 

M, 

dE (/ , 

E\ v 

{X d ,-X'd,)I dl 

+ E r 



dt 

T chi 

T do. 

T do, 



dE' ji 

E'd, 

{X ¥ -X' ¥ )I ¥ 





T qoi T qoi 


^ jdl _ Kjv + A:' | V m 

dt T m T m 


dVj, 

dt 


V, 


lii 


+ 


T T 

1 At 1 




\fi 


T T 

1 At 1 Fi 


T 


Ai 


jg t h , Ej± e 

dt T,, (T,f m 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
(2.11) 
(2.12) 

(2.13) 

(2.14) 


2. Algebraic Equations 


The algebraic equations consist of the stator algebraic equations and network 
equations. The stator algebraic equations directly follow from Fig. 2.1, applying 
Kirchhoff’s Voltage Law (KVL). A phasor diagram pertaining to universal rotating frame 
and local machine is given in Fig. 2.3. 
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Fig 2.3: Phasor diagram for universal rotating frame 

(a) Stator Algebraic Equations for i th machine, it is 

E\„ - V, sm(S, -0,)- R s ,I' h - = 0 

E ql - V l cos (S, -0,)- R Sl I qi - X'al a = 0 (2.15) 

(b) Network Equations 

Network equations corresponding to Fig. 2.4 are given as 
For Generator buses: 


Y,e‘°U a ~ jI v )e- Jls ‘-* ,3 > + P I AV,)+ jQu(V,) = £ v P k Y ik e > (2.16) 


k~\ 


i=l m 


For load Buses: 

R,AV,)+ JQ,AV,) = i V i V k Y, k e“< , ‘- B '- a » ) 


(2.17) 


k ~ 1 


i=m+l 
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Fig 2.4: Interconnection of synchronous machine dynamic circuit and the rest of the 
network 

In eq. (2.17) V^°' (Idrjlqi) e~ J (5, ~ k/2) =P Gi +jQa is the complex power “injected” into 
bus i due to generator. So Eq. (2.16) and (2.17) are only the real and reactive power 
balance equations at all the n-buses. Power balance equation at the generator buses shows 
the interaction of the state variables. 

Functionally the differential equations Eqs (2.8) - (2.14), together with stator 
algebraic equation. Eq.(2.15). and the network equations Eqs. (2.16)- (2.17) form a set of 
differential algebraic equations (DAE) of the form: 

x=f(x,y,u ) (2.18) 

0 = g(x,y) (2.19) 

With, x= [x/, X7, .... x„!]‘ as the state vector, where x, = [E q „ E S i: <w,-, Rf„ V Rl j, 

y~ fl'd-g, V, 6'] 1 as the set of algebraic variables where Id- q = [ Idi I q i Em E/mJ, and V= 

[V, V n ]\ 6= [6, 6 f j‘. 

As the network and stator transients are neglected, the current variables are of less 
important. For stability analysis, they are expressed in terms of the state and network 
variables. Either polar form or rectangular form can be used [21]. The network equations 
(2.16) and (2.17) are written in power balance form. They can also be expressed in current 
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balance form (. I=Y N V where Fa/ is N x N admittance matrix of the network). In transient 
stability analysis Fv is reduced from the standard load flow admittance matrix to take into 
account the machine internal reactance, and also the loads are represented as constant 
admittance type and included in reduced admittance matrix. This reduces the computational 
time and simplifies the analysis [21]. 

2.4 System Linearization 

Linearized analysis of multi machine power system is necessary for the study of both 
steady state and voltage stability. The DAE equations Eqs. (2.8)- (2.19) are linearized for 
analysis of low frequency oscillatory stability of the system. Two kinds of analysis are 
possible, viz 

(1) A multi machine linearized analysis that computes the eigen values and also finds 
those machines that contribute to a particular eigen value(participation factor analysis [ ]). 
Both local and inter area oscillations can be studied in such a framework. 

(2) A Single Machine Infinite Bus (SMIB) system case that investigates only local 
oscillations. 

Considering the participation of the generators in the oscillatory modes the site for 
PSS installations can be determined. The linearized analysis of the system is the basis for 
the designing of PSS for the machines [2], 

Considering Eqs (2.18) and (2.19), the vector y includes both the Ij. q and V vectors. 
Thus, Eq (2.18) is of dimension 7m, and Eq (2.19) is of dimension 2(n+m). To show 
explicitly the traditional load-flow equations and the other algebraic equations, y is 
partitioned as 

Y=[f d . q 6,V, v m \e 2 6 n V m+I VrJ ' 

= [y'a\y\] t (2.20) 

Here, the vector jy* corresponds to the load- flow variables and the vector y a correspond 

to the other algebraic variables. Bus 1 is taken as the slack bus and buses 2, m are the 

PV buses and my l, , n are PQ buses. The dimension ofx is 7m. Linearizing Eq(2.18) 

about a steady state operating point gives 
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Ax = AAx + BAy + EAu 
0 =CAx + DAy 


( 2 . 21 ) 

( 2 . 22 ) 


where, 


A = 


D = 


5f_ 

5x 
S g 


(x B ,y a ,u 0 ) 


8 y 


B = 


E = 


8 f 


J ( x „ , y „ , u „ ) 


8 y 
8 f 


C 


8 u 


J ( * „ , y „ , u „ ) 


5J_ 

8 x 


(*o»^o»«o) 

(2.23) 


The DAE equations can be expressed as 


Ax 



~A 

B 


Ax 

+ 

~E 

0 


_C 

D 




0 


(2.24) 


Matrix D is called the network algebraic jacobian [20], the system matrix (A sys ) is 
obtained after rearranging the algebraic variables [20]. 


Ax = A s y s A x + EAu 

(2.25) 

where, 


A. m = A - BD ~'C 

(2.26) 


The jacobian matrix gives the eigen values of the system for a small disturbance about 
the steady state operating point. The system is unstable if some of the eigen values lie in the 
right half of s-plane. For a weak system with increase in loading, the eigen values move 
towards the right half of s-plane [2], and after a particular critical value of loading the 
system becomes unstable. Also in some cases, the eigen value analysis shows that though 
the system is stable at steady state but for small disturbances the system starts oscillations. 
Oscillations in a system are undesirable as it affects the power transfer limit and also some 
cases sustained oscillation may lead to system instability. These small frequency 
oscillations are mainly because of fast acting voltage regulators and weak transmission 
lines (long lines). PSS is, thus, provided to damp out the oscillations in the system as 
described in chapter- 1. 
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2.5 Power System Stability and Power System Stabilizer 

Power System Stabilizer (PSS) is a supplementary control device to extend stability 
limits by modulating generator excitation to provide damping to the oscillations of 
synchronous machine rotors relative to one another. These oscillations of concern typically 
occur in the frequency range of approximately 0.2 to 2.5 Hz [19-21], and insufficient 
damping of these oscillations may limit the ability to transmit power. To provide damping, 
the stabilizer must produce a component of electrical torque on the rotor which is in phase 
with speed deviation [2]. For this PSS uses the as input the auxiliary signals which are 
derived from generator speed or electrical power output or field voltage etc [20]. The art 
and science of applying power system stabilizers have developed considerably, which 
involve the use of various tuning procedures and input signals, and learning to deal with 
practical problems such as noise and interaction with turbine generator shaft torsional 
modes of oscillation [2], The implementation details differ depending upon the stabilizer 
input signal employed. However, for any input signal the stabilizer must compensate for the 
gain and phase characteristics of the excitation system, the generator and the power system 
which collectively determine the transfer function from the stabilizer output to the 
component of electrical torque which can be modulated by excitation control [19]. This 
transfer function which can be denoted as GEP(s) is strongly influenced by voltage 
regulator gain, generator power level, and AC system strength. 



Fig. 2.5: Stabilizer with speed input-system block diagram 
The block diagram in Fig. 2.5 illustrates, in terms of a few basic small-signal transfer 
function blocks, the relationship between the applied torque on the turbine-generator shaft 
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and the resulting generator rotor speed, wc and rotor angular displacement §. The electrical 
torque may be considered to have two components, (a) that which is produced by the power 
system stabilizer solely by modulation of the generator flux, T ep , and (b) that which results 
from all other sources including shaft motion, T eo . The functional block relationship 
between speed and torque is shown for a stabilizer employing speed as an input signal. The 
contribution of torque due to the stabilizer path is given by: 

PSS a (s)* GEP (s) 

Acd g 

The transfer function GEP(s) represents the characteristics of the Generator, the 
Exciter, and the Power system. The variation of GEP(s) with exciter gain, generator 
loading, and ac system strength plays a dominant role in power system stabilizer tuning 
requirements and performance. The basic characteristics of this plant which are significant 
to stabilizer applications are as follows [2]: 

• The phase characteristics of GEP(s) are nearly identical to phase characteristics of 
closed loop voltage regulator. 

• The gain of GEP(s) increases with generator load. 

• The gain of GEP(s) increases as the ac system becomes stronger. This effect is 
amplified with high gain regulator. 

A properly tuned PSS must provide the necessary positive damping torque over a 
wide range of operating condition and disturbance to cancel the negative torque introduced 
by various sources [1], The overall excitation control system (Exciter, AVR, PSS and other 
Limiters) is designed so as to [4] 

• Maximize the damping of the local plant mode as well as the inter area mode 
oscillations without compromising the stability of other modes. 

• Enhance the system transient stability. 

• Minimize the consequences of excitation system malfunction due to component 
failures. 
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2.6 PSS Structure 

The block diagram of PSS, which is commonly used for analysis purposes, is given in Fig. 
2. 6. The PSS can be single input type or dual input type [10]. In case of dual input type, the 
blocks of PSS remain same for both of the inputs, only block parameters vary. The function 
of each block is given below. 



Fig. 2.6: PSS block diagram 


1. Washout Circuit 

The washout circuit is provided to eliminate steady state bias in the output of PSS, 
which will modify the generator terminal voltage. The PSS is expected to respond only to 
transient variations in the input signal and not the DC offsets in the signal. This is achieved 
by subtracting from it the low frequency component of the signal obtained by passing the 
signal through a low pass filter. 

The washout circuit acts essentially as a high pass filter and it must pass all 
frequencies that are of interest. If only the local modes are of interest, the time constant T w 
can be chosen in the range of 1 to 2 second. However if interarea modes are also to be 
damped, then T w can be chosen in the range of 10 to 20 second. There is noticeable 
improvement in first swing stability when T w is increased from 1.5 to 10. 

2. Dynamic Compensator 

The dynamic compensator used in practice is made up of multiple stages lead-lag 
compensators. A two stage lead-lag compensator has the following transfer function. 

T ( (j ) = Hill — ^ ^*1 ) ( ' + £ ^ 2 ) /q 

(1 + sT 3 )(1 + sT 4 ) ^ " 
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where Ks is the gain of PSS. The time constants 7/ to T 4 are chosen to provide a phase 
lead to the signal in the range of frequencies that are of interest (0.1 to 3 Hz). With static 
exciters, only one lead-lag stage may be necessary. In general, the dynamic compensator 
can be chosen with the following transfer function. 

K s N(s) 


T (s) = 


D(s) 


(2.28) 


where, 

N(s) = 1 +ajs+a 2 s 2 + +a n s n 

D(s) = l+bis+b 2 s 2 + +b p s p (2.29) 


The zeros of D(s ) should lie in the left half of s-plane. They can be complex or real. 
Some of the zeros of N(s) can lie in the right half of s-plane making it a non-minimum 
phase system. For design purposes, PSS transfer function is approximated to T(s), the 
transfer function of the dynamic compensator. The effect of washout filter can be neglected 
in the design, but must be considered in evaluating performance of PSS under various 
operating conditions. 


3. Limiter 


The output of the PSS must be limited to prevent the PSS action to encounter the 
action of AVR. For example, when load rejection takes place, the AVR acts to reduce the 
terminal voltage whereas PSS action calls for higher value of terminal voltage. It may even 
be desirable to block the PSS action in case of load rejection. The negative limit of PSS 
output is of importance during the back swing of rotor. The AVR action is required to 
maintain the voltage after angular separation has increased PSS action in the negative 
direction must be curtailed more than that in the positive direction. Ontario hydro uses a - 
0.05 p.u as the low limit and 0.1 to 0.3 as the higher limit [2], 

.4. Control Input 

The most commonly used and the obvious choice is deviation in the rotor angular 
speed. However, for practical implementations other signals such as bus frequency, 
electrical power, acceleration power and field voltage etc are also used [2]-[4], 
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2.7 Multi Machine PSS Tuning 


In case of PSS tuning considering SMIB model, the dynamic interaction between 
generators is not considered. This approach is sufficient when local mode of oscillation is 
of major concern. But when both local and inter area modes are considered, SMIB methods 
may not be sufficient. In usual practice PSS is tuned for local oscillations and when it is put 
into operation, the parameters are again adjusted depending on the oscillations with respect 
to other generators [21]. In case of multi machine PSS tuning methods, dynamic 
interactions between the generators are considered. This method is also called as 
coordinated calculation of PSS parameters. It is not necessary that all the machines will 
have the PSS; this method can also be applied when some of the generators are tuned for 
PSS depending upon the participation in oscillations. 

\2.7.1 Co-ordination of Power System Stabilizer Parameters 

This is called co-ordination method because the PSS parameters for all the machines 
are determined simultaneously. For a PSS with two lead-lag stages the PSS transfer 
function G pss (s) for a single input signal is given as 


G m (s) = 


Ar,(i + *r,)(l + ^) 
(\ + sT 2 )(l + sT A ) 


(2.30) 


where / and j are defined according to the signal used as input to the stabilizer, time 
constants T 2 and T 4 are assumed known, and the washout time constant T w will be 
considered in all cases to be equal to 10.0 second. The right hand side of (2.30) can be 
written as 

G pss (s) = d + —Jh±±JLQ ( 2.31) 

s + a 1 s + a 0 


where, 


a 1 


73+74 

t 2 t 4 


a 0 


1 KjTjTj 

T 2 T A 
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(2.32) 


P o 


K ‘-v -Pp)-A =-£rlTi+Tj -#(r 2 + r 4)] 


T 2 T 4 T 2 7V 


W 




An adequate state space representation for Gpss(s) can be obtained by using the 
observable canonical form as 


X c = A c x c + B c u c 


C r x ~ + D r u 


where 
Ac = 


0 

1 


a 0 

«1 


Be 


Po 

Pi 


,C c =[0 1 ],Dc = d 


(2.33) 

(2.34) 

(2.35) 


For the machine with a dual input PSS, Gpss(s) will be a 1 x 2 transfer function 
matrix of the form: 

Gpss(s) = [G,(s) G 2 (s)J (2.36) 

Where Gj(s) and G 2 (s) have the structure as given in equation (2.31). The time 
constants 7\ and Tj are considered to be the same for both signals. This means that matrices 
Ac and Cc will be the same. Matrices Be and Dc on the other hand, are augmented with a 
new column corresponding to the added input signal. 


2.7.2 Determination of Coordinated PSS Parameters 


In literature many methods are available for simultaneous tuning of PSSs. In this 
thesis two methods are used for this purpose. 

2.7.2. 1 Determination via Structurally Constrained Optimal Control [5], [6] 

PSS parameters can be obtained by using optimal controller design approach for a 
large power system. But these methods suffer from the disadvantage of needing the 
measurement and transmission of a large number of variables. This is a serious drawback 
concerning practical aspects. Again it can not be designed as a state feedback controller as 
all the states are not measurable especially in the case of power system. Constraints like 
output feedback (only measurable states are used) and decentralization (only local states are 
used) are of particular interest. 
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The PSS design problem consists of determining the parameters of the Gpss(s) transfer 
function so that eigen values of the compensated system are sufficiently stable. With T 2 and 
T 4 known, the only parameters which remain to be calculated are the PSS gain Kj and the 
lead-lag time parameters T, and 7}. 

Let the linear state space model of the power system without compensation is given as 
[5], [19]-[21]: 

x = Ax + Bu (2 37) 

y ~ C X (2.38) 

Compensating the system with the controller given by the Eq. (2.33) and (2.34) in 
feedback connections implies that u c = y and u=y c . The resulting composite system is then 
represented as: 


x =(A + B D c C)x + BC c x c 
x L = (A + BD l C)x + BC c x c 


(2.39) 

(2.40) 


Defining X a = [x 1 , x c r J 7 as the augmented state vector, Eqs (2.37) - (2.40) become 


x a — 


A + BD c C 
BC c C 


BC c 


(2.41) 


If in addition, we define a 3 x 1 augmented control vector it can be easily verified 
that Eq.2.23 can be rewritten as 

A " a = ( A a — B a G a C a )x a = A c \x a (2.42) 

where, 


Ga‘ 


A B C c 
0 A c 

- Dc T -Be 


B a = 


B 

0 


0 

I 


,C a = [C 0 ] 


(2.43) 


From the above equations, following informations can be inferred 
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y Matrices A a, Ba , and Ca are completely known as these consist terms T 2 and T 4 . 
The unknown parameters (J3 0 , (3\ and d, which depend on K h T, and 7}) appear only 
in matrix G a . 

y Eq. (2.24) exhibits the form of a constant-output feedback problem, where the 
(fictitious) process to be compensated is described by matrices A a, Ba, and Ca, and 

Ga is a constant output feedback matrix to be determined so that u a — -Gay a , where 
ya ~ C a x a . 

> Several methods can be used to determine the output feedback matrix G a , such as 

pole placement via output feedback or optimal output feedback, the later approach 
is preferred. 

y Some practical aspects of the multi machine problem impose further constraints on 
the control strategies. In fact, the conventional PSSs use signals derived from 
locally measurable quantities such as rotor speed, AC bus frequency or electric 
power. The use of cross feedback among the distinct machines, although 
theoretically desirable, poses practical difficulties. In a typical power system 
exhibiting dynamic stability problems, the power plants may be located hundreds of 
kilometers apart from each other. But the approach to process the information from 
different subsystems in order to tune PSS requires fairly sophisticated data 
transmission and coordination schemes. However, conventional PSS without using 
any cross-feedback can be tuned yielding good damping for all the modes. In order 
to achieve that, proper structural constraints like output feedback and 
decentralization are considered. 

With the system to be controlled given by Eqs.(2.37) -(2.38) , the objective is to 
determine the control strategy that minimizes the quadratic performance index: 

J ( x a , u a ) = ~df \ ( x a Q x a u a R u a ) d t (2.44) 

where the positive semidefinite matrix Q and the positive definite matrix R are state 
and input weighting matrices respectively. This is the so-called linear quadratic regulator 
(LQR) problem [ ] and its well known solution is given by the state feedback strategy: 
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u*=-K*x a (2.45) 

where, K = R 'Ba 1 P , and P is the solution of the Algebraic Riccati Equation (ARE) 

aI P + PA a - PB a R-'B T a P + Q = 0 (2,46) 

The above solution takes into account the whole system dynamics, as described, but it 
is again impractical concerning power system applications, since it prescribes feedback of 
al the state parameters to the entire input terminal. 

Applying the output feedback constraint, the control strategy can be expressed as: 

Wa G cl y u G a CaX a (2.47) 

That is, only the measured outputs are used to synthesize the system input. 
Comparing this with Eq. (2.45) by imposing the output feedback constraint on the 
conventional feedback strategy is tantamount the following equation for Ga : 

£*= G a C a (2.48) 

t* t* ♦ t* nr 

Since above equation implies that C a G a =K , matrix G a can be obtained by using 
the concept of matrix pseudo inverse, applied to matrix Cj . It can, thus, be easily shown 
that the output feedback gain matrix is given by 

G a = KCj (CaCjy 1 (2.49) 

As K matrix is not a full matrix, the pseudo inverse is taken, for which, this method is 
called “Pseudo decentralization method” 

Now applying the decentralization constraint that is only the local states are fed back 
to the local inputs. It is easy to see that the decentralization constraint imposes a block 
diagonal structure on the state feedback matrix K. Therefore, the decentralized can be 
expressed as: 

K-Block diagonal {K1 K 2 K 3 K N } =0 (2.50) 

The generalized ARE considering both the constraints (output feedback and 
decentralization) is given as 
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a£ P + PA a - PB a R ] b£p + Q+L T RL = 0 (2.51) 

where L matrix takes into account the two constraints and determined using in 
iterative process. The multiplying factors selection is quite arbitrary so also selection of 
weighting matrices Q and R. Initially P is calculated by solving the generalized ARE and 
using this value, K is calculated. From K matrix, Ga is calculated. Knowing the Ga matrix 
the block diagonal elements give the PSS settings, K h Ti and 7} 

2.1.22 Determination by Strip Eigen Value Assignment Method [7], [8] 

The problem with the pseudo decentralization method is that the weighting matrices 
Q and R are selected arbitrarily. In some literature [6], they are taken as identity matrices. 
In this case, the PSS gives good damping for a particular operating condition only. But for 
some other loading conditions the PSS puts adverse effect on the eigen values, some 
become negative also. Proper weighting matrix selection takes a large time which is not 
suitable for real time application. Some literatures [6] have taken R as an identity matrix 
and values of diagonal elements of Q depending upon the damping of that particular mode. 
Here we have tried with giving weight to all the modes arbitrarily which leads to an optimal 
setting for many operating conditions. 

In the strip eigen value assignment method there is no need to define the Q matrix at 
any stage. A linear quadratic controller is designed such that the optimal closed-loop 
system has eigen values lie within a vertical strip in complex s-plane [8]. The desired 
position of eigen values is obtained without any convergence problem. 

Let (A, B) be the pair of the open-loop system matrices in the system Eq. (2.37) and 
h(>0) represents the prescribed degree of relative stability. Then the closed loop system 
matrix Ac-A-BR'' B 7 P has all the eigen values lying on the left side of the -h vertical line as 
shown in Fig. 2.7 where matrix P is the solution of ARE 

(A + hlfP+P (A +hln) - PBR' 1 B T P+Q~O n (2.52) 

With (>=0„, the unstable modes of A+hI n are shifted to their mirror image positions 
with respect to -h vertical line, which are the eigenvalues of the closed loop system matrix 

Ac- 
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Fig. 2.7: (a) single vertical line (b) vertical strip 

Assume that hi and h 2 are two positive real values to determine an open vertical strip 
of [-hi, -h\\ on the negative real axis as shown in Fig. 2.7 and given as an n x n matrix Ai = 
A +h/I. The control law is changed to be 

u (0 = -pFx (0 (2.53) 

With the feedback gain F=R' I B I P. The matrix P is the solution of the following 
modified ARE: 

A i T P+ PA 1 -PBR ' 1 B r P~0„ (2.54) 


P 


The gain p is selected by 
h 2 - h\ 


0.5 + 


2tr(Aj ) 


0.5 + 


h 2 ~ h\ 
tr(B.F ) 


(2.55) 


where trace (Ai) are the eigen values of A/ in the right half of the complex s-plane. 
The final closed loop system consists of a set of eigen values which lie inside the vertical 
strip of [-h 2 , -hi). Flere we can assume the matrix R to be a unity matrix. Thus solving the 
ARE does not require Q matrix. 


2.8 Results and Discussions 


The pseudo decentralization method and the strip eigen value placement method have 
been implemented on WSCC 3-machine, 9-bus system and New England 39-bus system. 
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Though, the objective is to tune PSS for NPC system, the standard systems were tested 
for the verification of proposed methods. The line diagram and the system data are given 
in appendix[l]. MATLAB® is used for all the analysis purposes 

Initially the base case parameters were found out for the given operating and loading 
conditions. After linearization of the DAE equations around the base case points the eigen 
values were calculated using the code developed in MATLAB®. Other parameters such as 
damping ratio, frequency of damped and undamped oscillations, and the participation 
factors were calculated. Participation factors decide the location of PSS. Then PSS was 
considered and the new eigen values were calculated. The zero eigen values appearing in 
the analysis are called as redundant eigen values. This is because there is no infinite bus, so 
the angles are relative angles and the machine inherent damping is taken as zero [19], [21]. 

2.8.1 WSCC 9-bus System 

WSCC system consists of three generators. Bus 1 is considered as infinite or swing 
bus for all the calculations. The generators are equipped with rotating type exciters. The 
rate stabilizer, AVR amplifier and exciter were considered in detail with the exponential 
saturation function [19]. The participation factor analysis shows that generator two and 
three are participating maximum in the critical mode shown in Table 2.1. The position of 
eigen values in s-plane are shown in Fig. 2.1. So PSS were placed at generators 2 and 3. The 
eigen value after inclusion of PSS shows improvement in system damping. Table 2.2 shows 
the eigen values after placing PSS at generators 2 and 3 using pseudo decentralization 
method. The eigen values with strip pole placement method after putting PSS at generators 
2 and 3 is given in Table2.3.As the objective is only to test the proposed method non-linear 
simulation was not done for the WSCC system. Both Pseudo decentralization and strip 
eigen value assignment methods are tested. Figs 2.10 and 2.11 shows the eigen values 
computed using strip eigen value placement method and pseudo decentralization method. It 
can be observed that with pseudo decentralization method the minimum damping ratio is 
0.0531 and with strip eigen value method 0.0985. Table 2.4 shows the parameters of PSS 
obtained by using pseudo decentralization method with speed deviation as input to PSS 
(single input PSS). 
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Table 2.1 : Eigen values of 3-machine 9-bus system 


Eigen values 

Damping ratio 

Frequency of oscillation in Hz 

0.0000 

0.0000 

0 

0.0000 

0.0000 

0 

-0.4488 ±0.5773i 

0.6137 

0.1164 

-0.4326 ± 0.7762i 

0.4868 

0.1414 

-0.4424 ± 1.2030i 

0.3451 

0.2040 

-3.2258 

1.0000 

0.5134 

-3.5139 

1.0000 

0.5593 

-4.9221 

1.0000 

! 

0.7834 

-0.2084 ± 7.5727i 

0.0275 

1.2057 

-5.0877 ± 7.7705 i 

0.5478 

1.4782 

-5.1683 ± 7.87361 

0.5478 

1.4990 

-5.1839 ± 7.9161 i 

0.5478 

1.5060 

-0.8377 ±11.49031 

0.0727 

1.8336 


Eigen Value location 



Fig. 2.8: Eigen values for 3 machine 9 bus system 





Table 2.2: Eigen Values with PSS at two machines (Pseudo decentralization method) 


Eigen Values 

Damping Ratio 

Frequency of Oscillation in Hz 

0.0000 

0 

0.0000 

-0.3204 ± 0.0628i 

0.9813 

0.0520 

- 0.5000 

1.0000 

0.0796 

- 0.5000 

1.0000 

0.0796 

-0.4274 ± 0.656 li 

0.5458 

0.1246 

-0.4441 ±0.8028i 

0.4841 

0.1460 

- 1.0000 

1.000 

0.1592 

- 1.0000 

. 1.000 

0.1592 

-1.0772 

1.000 

0.1714 

-0.2689 ± 1.20731 

0.2174 

0.1969 

' -3.2258 

1.0000 

0.5134 

-3.5 1 28 

1.0000 

0.5591 

-4.0187 

1.0000 

0.7828 

-0.4029 ± 7.57421 

0.0531 

1.2072 

-4.9446 ± 7.6568i 

0.5425 

1.4506 

-5.1376 ±7.90851 

0.5448 

1.5010 

-5.2122 ±7,86771 

0.5523 

1.5020 

-0.865 1 ±11 .50791 

0.0750 

1.8367 




Table2.3 Eigen value for 3- machine 9-bus system with strip eigen value placement method 


Eigen V alue 

Damping ratio 

Frequency of Oscillation in 

Hz 

-0.4845 ± 0.5454i 

0.6641 

0.1161 

-0.4548 ± 0.7 154 i 

0.5365 

0.1349 

-0.4628 ± 1 . 1 545i 

0.3721 

0.1980 

-1.7056 ± 0.6927i 

0.9265 

0.2930 

-3.2258 

1.0000 

0.5134 

-3.4975 

1.0000 

0.5567 

-4.9032 

1.0000 

0.7804 

-1.5568 ± 7. 63 98 i 

0.1997 

1.2409 

-5.2324 ± 7.8321 i 

0.5555 

1.4991 

-5.4004 ± 7.91 79 i 

0.5635 

1.5254 

-5.5650 ± 7.949 li 

0.5735 

1.5444 

-1.1365 ±ll.4838i 

0.0985 

1.8366 


Eigen value location 



Fig. 2. 1 0: Eigen Value of 3 machine 9-bus system with strip eigen value placement 


Table 2.4 PSS parameters for 3 machine 9-bus system with pseudo decentralization 


Machine No. 

K 

T. 

t 2 

t 3 

t 4 

2 

2.213 

0.1525 

0.0525 

0.05913 

0.05 

3 

1.49 

0.4765 

0.05 

0.0210 

0.05 
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2.8.2 New England 39-bus system 

I he New England system consists of 10 machines and again the generators are 
equipped with DC exciters. Ihe machine at bus 1 is considered are infinite generator for all 
the analysis. Ihe open loop analysis (Table 2.5) shows that, there are three critical modes 
(modes with damping ratio < 0.05). Fig. 2.12 shows the eigen values for open loop case. 
Though the generators 3 and 7 do not participated in the oscillations corresponding to these 
modes, PSS was assumed in all the 9 generators except the 1 st generator. The strip eigen 
value placement method can not be applied to this system as it has got some limitations for 
bigger systems. A dual input PSS was considered. The parameters after inclusion of PSS 
are given in Table 2.6. It shows marked improvement in the eigen values. Fig.2.12 shows 
the eigen value location after inclusion of PSS. The minimum damping ratio after inclusion 
of PSS is 0.0528. The improvement from previously developed method is that Q matrix 
was assigned values depending on the participation of the generators in the oscillations 
[19]. 


Eigen value plot 



Fig. 2.11: Eigen values for 1 0 machine 39-bus New England system 
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Table 2.5 Eigen values of 10 machine 39 bus system 


Eigen Value 

( 1 .0e+002 * ) 

Damping 

ratio 

Frequency 

(Hz) 

Eigen Value 

( 1 .0e+002 * ) 

Damping 

ratio 

Frequency 

(Hz) 

0.0000 

0 

o.oooo 

-0.0449 

1.0000 

0.7150 

-0.0000 

0 

0.0006 

-6.0507 

1.0000 

0.8063 

-0.0040 

1.0000 

0.0636 

-0.0020 ± 0.0548i 

0.0367 

0.8728 

-0.0027 ± 0.004 li 

0.5481 

0.0785 

-0.0046 ±0.0559i 

0.0814 

0.8922 

-0.0057 

1.0000 

0.0914 

-0.0026 ± 0.05991 

0.0441 

0.9538 

-0.0057 ± 0.00 18i 

0.9556 

0.0949 

-0.0629 

1.0000 

1.0007 

-0.0047 ± 0.0037i 

0.7809 

0.0949 

-0.0041 ± 0.0667i 

0.0610 

1.0630 

-0.0047 ± 0.0049i 

0.6895 

6.1080 

-0.0708 

1.0000 

1.1264 

-0.0071 ± 0.0053i 

6.8041 

6.1415 

-0.0035 ± 0.0720i 

0.0479 

1.1472 

-0.0083 ± 0.00 5 9 i 

0.8145 

0.1622 

-0.0745 

1.0000 

1.1855 

-0.0060 0.00891 

0.5596 

0.1717 

-0.0048 ± 0.0778i 

0.0619 

1.2410 

-0.0109 ± 0.00061 

0.9985 

6.1734 

-0.0068 ± 0.0847i 

0.0795 

1.3523 

-0.0138 

1.0000 

0.2204 

-0.0057 ± 0.0857i 

0.0662 

1.3663 

-0.0144 

1.0000 

0.2295 

-0.1558 

1.0000 

2.4789 

-0.0156 

1.0000 

0.2485 

-0.1573 

1.0000 

2.5043 

-0.0089 ± 0.0 1 29i 

0.5688 

0.2501 

-0.1576 

1.0000 

2.5076 

-0.0166 

1.0000 

0.2639 

-0.1759 

1.0000 

2.7999 

-0.0178 ±0.0048i 

0.9656 

0.2933 

-0.4827 

1.0000 

7.6822 

-0.0205 ± 0.0017i 

0.9965 

0.3268 

-0.4845 

1.0000 

7.7104 

-0.0066 ± 0.01981 

0.3151 

0.3322 

-0.4915 

1.0000 

7.8221 

-0.0259 

1.0000 

0.4117 

-0.4926 

1.0000 

7.8395 

-0.0271 ± 0.0083i 

0.9559 

0.4511 

-0.4931 

1.0000 

7.8471 

-0.0019 ±0.0369i 

0.0522 

0.5886 

-0.4990 

1.0000 

7.9418 

-0.0443 

1.0000 

0.7049 

-1.1688 

1.0000 

18.6023 
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Table 2.6 Eigen values ol New England 3 9 -bus system with Pseudo decentralization 


Eigen Value 

Damping 

Frequency 

Eigenvalue 

Damping 

Frequency 

(1.0e+002*) 

Ratio 

in Hz 

(1.0e+002*) 

Ratio 

in Hz 

0.0000 

0 

0 

-0.0446 

l.oooo 

0.7101 

-0.0023 

1 .0000 

0.0361 

-0.0498 

1.0000 

0.7919 

-0.0040 

1.0000 

0.0631 

-0.0044 ± 0.0555i 

0.0793 

0.8865 

-0.0026 ± 0.004 li 

0.5364 

0.0782 

-0.0058 ± 0.05691 

0.1018 

0.9097 

-0.0058 

1.0000 

0.0921 

-0.0035 ± 0.0601 i 

0.0579 

0.9575 

-0.0054 ± 0.0025i 

0.9083 

0.0951 

-0.0636 

1.0000 

1.0124 

-0.0046 ±0.0040i 

0.7551 

0.0970 

-0.0057 ±0.0675i 

0.0837 

1.0779 

-0.0047 ± 0.0052i 

0.6680 

0.1122 

-0.0711 

1.0000 

1.1320 

-0.0060 ± 0.0053 i 

0.7920 

0.1387 

-0.0038 ± 0.0722i 

0.0528 

1.1500 

-0.0062 ± 0.00 86i 

0.5858 

0.1691 

-0.0743 

1.0000 

1.1828 

-0.0108 

1 .0000 

0.1712 

-0.0049 ± 0.0779i 

0.0626 

1.2416 

-0.0074 ± 0.008 li 

0.6734 

0.1749 

-0.0079 ± 0.0847i 

0.0923 

1.3546 

-0.0117 

1.0000 

0.1861 

-0.0073 ±0.0866i 

0.0836 

1.3828 

-0.0123 

1 .0000 

0.1960 

-0.1573 ± O.OOOli 

1.0000 

2.5035 

-0.0142 ± 0.0005 i 

0.9995 

0.2259 

-0.1580 

1.0000 

2.5146 

-0.0157 ±0.0002i 

0.9999 

0.2504 

-0.1758 

1.0000 

2.7985 

-0.0049 ± 0.01 62 i 

0.2902 

0.2699 

-0.2000 

1.0000 

3.1831 

-0.0193 

1.0000 

0.3070 

-0.2000 ± O.OOOOi 

1.0000 

3.1831 

-0.0199 ± 0.0022i 

0.9937 

0.3193 

-0.2000 ± O.OOOOi 

1.0000 

3.1831 

-0.0207 ± 0.008 li 

0.9310 

0.3543 

-0.2000 

1.0000 

3.1831 

-0.0250 

1 .0000 

0.3972 

-0.2000 

1.0000 

3.1831 

-0.0028 ± 0.028 li 

0.1001 

0.4491 

-0.2000 

1.0000 

3.1832 

-0.0400 

1 .0000 


-0.2000 

1.0000 

3.1835 

-0.0400 ± O.OOOOi 

1 .0000 


-0.2001 

1.0000 

3.1847 

-0.0400 

1.0000 


-0.4840 

1.0000 

7.7024 

-0.0400 ± O.OOOOi 

1.0000 


-0.4856 

1.0000 

7.7293 

-0.0400 ± O.OOOOi 

1.0000 





-0.0401 

1.0000 

0.6383 

-0.4926 

1.0000 

7.8400 

-0.0405 

1.0000 

0.6449 

-0.4931 

1.0000 

7.8478 

-0.0083 ± 0.04 14i 

0.1970 

0.6724 

-0.4990 

1.0000 

7.9421 

-0.0423 

1.0000 

0.6733 

-1.1688 

1.0000 

18.6025 
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Eigen value plot 



Fig. 2. 12: Eigen Value for 10 m/c 39 bus system with Pseudo Decentralization 
2.9 Conclusion 

The methods described in this chapter can be invariably used for PSS tuning in case 
of multi machine environment. But still the analysis is not complete. Before doing multi 
machine PSS tuning we must first perform the SMIB analysis, because a lot of literatures 
are available in this respect and till date PSS tuned through SMIB are operating well. But 
the truth is that, we should consider the dynamics among all the generators, this will reduce 
the burden of doing trail and error methods in the field while practically incorporating the 
PSS. Again strip eigen value method and pseudo decentralization method are capable of 
shifting the under damped eigen values (Eigen values with £ <0.05) to the safe damping 
region. But strip eigen value placement has got inherent disadvantage of application to 
small system only. The selection of the weighting matrices is still a hit and trial method. In 
some cases we take them as identity matrices [6], and in some cases we are giving 
weightage depending on its participation in system oscillating modes. Depending on the 
operating conditions also these matrices has to be changed which is practically very 
difficult to achieve. This process is also an arbitrary method. Some robust intelligent PSS 
has to be found out, which will be simpler and easy to implement. 
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Chapter 3 


PSS Tuning for Kaiga Nuclear Power Plant 

3.1 Introduction 

Kaiga nuclear power plant has static busfed type of exciter, which gets power from an 
auxiliary transformer. The basic function of AVR is to control the firing angle of the 
exciter. So AVR output is fed to the exciter gate circuit. From the AVR output to the main 
held of generator a linear circuit is considered [19] i.e. the rectifier dynamics are neglected. 
For fast operation during disturbance the AVR regulator is provided with PID 
characteristics which do not need a rate feedback [19]. Also it is provided with high ceiling 
voltages for field forcing. AVR time constant and gain can be varied by selector switches in 
the plant. Fig. 3.1 shows the block diagram of AVR. The current settings and range of 
values are given in Appendix-C. 



Vo*(l+sT,) 

1+sWVp) 



Fig. 3.1 : AVR block diagram 
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3.2 Exciter and AVR Dynamic Equations for Kaiga NPP 

Kaiga AVR has no feedback stabilization (rate feed back) and the PID amplifier has 
equal number of poles and zeros, this is reduced in the form of observable canonical form 
[6]. There is no saturation as the exciter is static type. The limits are not considered for 
linear modeling. Fig.3.2 shows the simplified AVR block diagram. The initial conditions 
for steady state load flow were calculated from the dynamic equations. 

The method is similar to that discussed in the section 2.6.1. From the Fig. 3.1, the 
output y c is the same as Vr (fig. 2. 2). Error voltage ( Vter- Vref) is the control input w Cj 


Uc V r ef “Vter ^ 

V 0 *Q.+sT a ) 


(1 +sT b ) 

yc=v R _ 

— — p. 

(l+sT a V 0 /V P ) 


(1 +sT b V M !V p ) 

w 


Fig. 3.2: Simplified AVR block Diagram 


The PID transfer function is given as 
K s (1 + sTj )(1 + sT i ) 

s } = (1 + 7t~ 2 ) (i + 7u~) 


(3-1) 


The state space model in observable canonical form can be written as 


x, 

*2 



1 

0 

1 

o 

L 


~X\ ‘ 


'fio~ 





4- 



1 - a\ 


_ x 2 _ 




[0 ij 


*1 

*2 


+ d u c 


(3-2) 

(3-3) 


a.o, at, fio, 0i, and d in terms of K s , 7), 7}, 73, and T4 are given in Eq (2.32). y c in terms 
of the actual parameters are given as, y c = Vr, u c = V re f- V ter . The initial conditions can 
be calculated by equating derivative terms to zero at steady state [20]. These equations are 
linearized and included for further analysis. 
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3.2.1 PSS Block Diagram at Kaiga NPP 

Kaiga NPP has got a dual input PSS as shown in Fig.3.3. Electrical power deviation 
and speed deviation or frequency deviation are the two inputs. The power and frequency 
transducers sense the actual values. After weighting and washout stage, deviation in 
electrical power and deviation in frequency signals are obtained. Depending on the power 
deviation, gain K 3 is decided. In this case, for power deviation less than 15% of rated value, 
fa varies as a quadratic function of AP, which implies that PSS is not fully utilized. For AP 
greater than 15% of rated value, K 3 is 1 that is PSS is fully utilized. For tuning the PSS, 
lead-lag block are not necessary to be tuned. Properly weighted (multiplied with K 3 ) 
frequency or speed deviation signal is added to negative of the weighted power deviation 
signal gives final PSS output. For proper phase compensation so that final torque due to 
PSS is in phase with the speed deviation signal, K/ and K 2 are to be determined. The range 
of the parameters is given in Appendix-C 



Fig. 3.3: PSS block diagram of Kaiga NPP 


3.3 SMIB Analysis and PSS Tuning 

The SMIB approach is used for tuning of Kaiga PSS. The reason behind this is, 
nearby hydro units are facing no problem when there is oscillation at Kaiga, so they do not 
have immediate plans to go for PSS tuning. The second reason can be attributed to the non- 
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availability of exact data of those generators. Even though there are different multi machine 
PSS tuning methods, SMIB method is still preferred as the other developed methods are 
very difficult to implement [21]. 

3.3.1 Phase Plot 

The single machine equivalent was obtained from the existing data for the two 
machines. The voltage and MVA base were changed. The two nearby generators were 
assumed as constant impedance loads injecting power to the system [21] as opposed to the 
normal loads. One infinite bus was assumed at Guttur connected to Kaiga over the double 
circuit line of 165 km length. After static load flow analysis, the loads were changed to 
equivalent admittance knowing the voltage magnitude at their terminals. Then Thevinin 
equivalent was calculated at the terminals of Kaiga, transfer reactance was calculated 
between the single machine equivalent and assumed infinite system. The constants K1-K6 

[1] , [10] were calculated. The transfer function for Generator, Exciter and Power system, 
GEP(s) [2] was found out. This is nothing but the transfer function from the PSS input 
point to the output of exciter as in the Heffron-Phillips model [19]. Phase plot which is 
given in Fig 3.5 and Gain plot which is given in Fig.3.6 were obtained for the system at 
different frequencies. This gives the amount of phase compensation required by the PSS 
output so that the final torque produced due to this will be in phase with the speed deviation 

[ 2 ] , 


phase plot 



Fig. 3. 4: Phase plot for kaiga NPP 
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Fig. 3.5: Gain Plot for KaigaNPP 
3.3.2 Eigen Value Analysis 

The linearized model as developed in 2.4 is replaced with Kaiga AVR and exciter. 
T he Inter Connecting Transformer (ICT) is considered as a part of the network for eigen 
value analysis. The two nearby hydro generators are considered as two loads (of 90 MW 
each) for this analysis. As they are small in size and do not oscillate when Kaiga generators 
arc oscillating, the SMIB analysis gives enough idea about the system behaviour for small 
disturbances. The base case load flow is run and the system is linearized. The eigen value 
analysis is performed and given in Table 3.1. The minimum damping ratio of the system is 
found to be 0.0102. Fig. 3.6 shows the eigen values in s-plane. Two critical modes are 
existing in the open loop case (without PSS) which proves the insufficient damping for the 
existing system. Participation factor analysis shows that Kaiga generators are participating 
maximum in the oscillatory modes corresponding to speed and load angle deviation (Aco 
and A5). The frequency of oscillation comes out to be in the range of local mode [19]. So it 
proves that Kaiga generators are oscillating with respect to a strong system connected over 
a long transmission line. No intraplant mode of oscillation is observed. The oscillation 
pattern changes for different loading and AVR setting conditions. So the cause of 
oscillation can be attributed to the improper setting of the AVR as after inclusion of ICT 
also there is no major improvement in the eigen values. Normally inclusion of ICT solves 
the problem of generator connected to rest of the system with weak line (Xe >0.5 PU, 
including transformer reactance). The zero eigen values appearing in the analysis are 
redundant eigen values as discussed in section 2.7. 
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3.3.3 PSS Tuning 


The PSS provided by BHEL has no lead-lag block as the conventional PSS has. So 
the methods used lor tuning PSS are not used here. The tuning is done completely through 
MATLAB 1 ' 1 ’ SIMIJLINK which is explained in the next section. The PSS has options for 
either frequency or speed input. Here we have used the direct speed input. Frequency input 
PSS is for weak to medium transmission systems [2]. As the phase plot shows that there is 
considerable phase lag in the range of frequency of interest, only local mode is addressed. 
The eigen value analysis confirms the oscillations at local mode and poor d am ping of the 
oscillations. As there is no inter area mode of oscillations SMIB analysis is sufficient even 
though this aspect also needs attention [21]. Ki and K 2 values are the important 
parameters, which can only be varied. Other parameters such as the transducer time 
constant, washout filter time constant and final gain are fixed values. Appendix [3]. After 
inclusion of PSS the system eigen values are given in Table 3.2. From Table 3.2 the 
minimum system damping is 0.0729. This shows the improvement in the system damping. 
Fig. 3.7 shows the system eigen values after PSS is included in the linearized analysis. 


Table3.1 Eigen value analysis without PSS 


Eigen Value( l.0e+002 *) 

Damping ratio 

Frequency in Hz 

0.0000 

0 


- 0.0000 

1.0000 


-0.0054 ± 0.0067i 

0.6308 

0.1371 

-0.0119 

1.0000 

0.1895 

-0.0108 ± 0.0059i 

0.8766 

0.1952 

-0.0264 

1.0000 

0.4202 

-0.0408 

1.0000 

0.6488 

-0.0479 

1.0000 

0.7617 

-0.0006 ± 0.0498i 

0.0102 

0.7920 

-0.0623 

1.0000 

0.9923 

-0.0045 ± 0.06741 

0.0666 

1.0752 

-03054 

1.0000 

4.8610 

-0.3139 

1.0000 

4.9966 

-0.3165 

1.0000 

5.0369 

-8.3033 

1.0000 

132.1519 

-8.3223 
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Fig. 3.6: Eigen value without PSS 


Table 3.2 Eigen Value Analysis with PSS on 1 st Machine at Kaiga NPP 


Eigen Value( i.0cK)()2 *) 

Damping ratio 

Frequency in Hz 

0.0000 

0 

0.0000 

-0.000 

1.0000 

0.0000 

-0.0017 

1.0000 

0.0266 

-0.0055 ± 0.00691 

0.6234 

0.1412 

-0.0096 

1.0000 

0.1533 

-0.0139 

1 .0000 

0.2205 

-0.0245 

1 .0000 

0.3901 

-0,0249 

1.0000 

0.3960 

-0.0492 

1.0000 

0.7827 

-0.0041 ± 0.0560i 

0.0729 

0.8943 

-0.0637 

1.0000 

1.0140 

-0.0356 ±0.20721 

0.1694 

3.3456 

-0.3081 

■ 

4.9034 

-0.3164 

1.0000 

5.0361 

-0.3232 

1 .0000 

5.1440 

-8.2554 

1.0000 

131.3884 

-8.3124 

1.0000 

132.2964 

-8.3281 

1.0000 

132.5456 
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Fig 3.7: Eigen value locations with PSS on generator-1 
3.4 Simulink Models for NPC System 

Non-linear Simulation is required for the transient stability analysis of the system. For 
this, Simulink 0 " model is used. It uses the generator model as described in section 2.2 with 
transient saliency neglected (i.e. X q = X d). Other detail system data are given in Appendix- 
3. Swing bus is designed with a generator having very high Inertia constant and large 
capacity so that rotor dynamics can be neglected. 

Simulink model of Fig. 3.8 shows the reduced model of NPC system, where the 
nearby hydro generators are considered as loads of 180MW. A large value of load is also 
assumed at the swing bus end. The AVR, exciter and PSS simulink models (Fig. 3.9, Fig. 
3.10, and Fig. 3.11) used are similar to those supplied by the manufacturer. Fig 3.12 shows 
the detail NPC system model. The two nearby generators are considered and one of them is 
equipped with IEEE-1 exciter. The NPC system model is further expanded by considering 
16 other generators to make the system more practical (Fig. 3.18). Some of these generators 
are equipped with exciters and others have constant field voltage. Internal subsystem 
diagram for one machine is shown in Fig. 3.14. Kaiga generators are equipped with IEEE-1 
exciter which is shown in Fig. 3.15, for a comparative study purpose. 
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Fig. 3.8: Simulink blocks of reduced NPC system 
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Fig.3.9: Exciter and AVR signal collection diagram 



Fig. 3.10: Exciter and AVR simulink model 
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Fig. 3.12: Detail model ofNPC system 
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Fig.3.14: Kaiga system with IEEE-1 exciter 
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3.5 Simulation and Results 


Using Giaphical Usci Interface (GUI) block, the initial steady states and load flow is 
perlormed. I his is necessary to start the simulation from steady state. A line is tripped at 30 
sec lor all the simulations. Kaiga-1 and Kaiga-2 are rated at 200MW and 215 MW 
respectively. Following cases were studied. 

(a) NPC reduced system without PSS 

(b) NPC reduced system with PSS at generator-1 

(c) Detailed NPC system without PSS 

(d) Detailed NPC system with PSS at generator-1 

(e) Reduced NPC system with 16 other generators 

(1) Reduced NPC system with 1 6 other generators and PSS at Kaiga-1 

3.5.1 NPC Reduced System without PSS 

The time constant ratio Tj /'IN is set to 0.5 and the simulation is performed with a fault 
initiated at 30 second. The load angle and power variations are shown in Fig. 3.15 and Fig. 
3.16, respectively. The speed deviation between the both Kaiga generators is plotted in Fig. 
3.17 whereas speed deviation between Kaiga and slack bus (bus-1) is shown in Fig. 3.18. 
The swing generator load angle and Kaiga bus voltage are also oscillating which can bee 
seen from Fig. 3.19 and Fig. 3.20 respectively. The frequency of oscillation is observed as 

1 .2 Hz which is quiet similar to the result of eigen value analysis. By the changing the Tj 
/T 2 ratio to 0.25, the system oscillations are getting damped after few cycles which can be 
seen from the Figs. 3.21 to 3.24. From Figs. 3.21 to 3.24 it can be observed that oscillation 
is damped out after more than 1 5 second. 

3.5.2 NPC Reduced System witli PSS at Generator-1 

PSS signal is given to the AVR summing point of gcn-1, the time constant ratio Tj /T 2 
is set to 0.5, and PSS settings K|=0.4, K2=0.5. Simulation is performed with a fault initiated 
at 30 second. The load angle and power variations are shown in Fig. 3.25 and Fig. 3.30 
respectively. The speed deviation between the both Kaiga generators is plotted in Fig. 3.26 
whereas speed deviation between Kaiga and slack bus (bus-1) is shown in Fig. 3.27. The 
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swing generator and Kaiga bus voltage are shown in Fig. 3.28 and Fig. 3.29 respectively. It 
can be observed that there is very little oscillation in system parameters. 

3.5.3 Detailed NPC System without PSS 

The two nearby hydro generators were connected to the Kiaga generators. The exciter 
parameters of Kadra hydro are assumed suitable values [15]. The time constant ratio Tj /T 2 
at Kaiga- 1 is set to 0.5 and the simulation is performed with a fault initiated at 30 second. 
The load angle of Kaiga-1 and Kudasalli hydro are shown in Fig. 3.31 and Fig. 3.33 
respectively. The speed deviation between the both Kaiga generators is plotted in Fig. 3.34 
whereas speed deviation between Kaiga and Kadra is shown in Fig. 3.32. The frequency of 
oscillation is observed as around 1.2 Hz which is quiet similar to the result of eigen value 
analysis. . From Figs. 3.31 to 3. 34 it can be observed that oscillation is not damped out. 

3.5.4 Detailed NPC System with PSS at Generator- 1 

PSS signal is given to the AVR summing point of gen-1, the time constant ratio Tj /T 2 
is set to 0.5, and PSS settings Ki=0.5, K 2 =0.8. The simulation is performed with a fault 
initiated at 30 second. The load angle of Kaiga-1 and Kadra hydro are shown Fig. 3.35 and 
Fig. 3.37 respectively. 'The speed deviation between the both Kaiga generators is plotted in 
Fig. 3.36 whereas speed deviation between Kaiga and Kadra is shown in Fig. 3.38. The 
PSS output and Kaiga-1 bus voltage are shown in Fig. 3 .39 and Fig. 3.40 respectively. It 
can be observed that there is very little oscillation in system parameters. 

3.5.5 Reduced NPC System with 16 Other Generators 

To make the system more practical 16 other generators are connected at the swing bus 
end. Some of these generators have exciter on them. The time constant ratio Tj /T 2 is set to 
0.5 and the simulation is performed. Following a three phase fault and subsequent opening 
of line connecting Kaiga and the nearby load, Kaiga-1 and 2 start oscillations. The load 
angle variations of generator- 1, generator at 200km from Kaiga, and swing generator are 
shown in Fig. 3.41, Fig. 3.42 and Fig. 3.45 respectively. The speed deviation between the 
both Kaiga generators is plotted in Fig. 3.43 whereas speed deviation between Kaiga and 
slack bus (bus-1) is shown in Fig. 3.46.Bus voltage at Kaiga-1 is also oscillating which can 
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bee seen from Fig. 3.44. From Figs. 3.41 to 3.46 it can be observed that oscillation is 
damped out after more than 25 second. 

3.5.6 Reduced NPC System with PSS at Kaiga-1 

PSS signal is given to the AVR summing point of gen-1, the time constant ratio Tj /T 2 
is set to 0.5, and PSS settings K|— 0.2, K 2 — 0.3. The simulation is performed with a fault 
initiated at 30 second. The load angle of Kaiga-1, generator at 200km from Kaiga and 
swing generator are shown Fig. 3.47, Fig. 3.48 and Fig. 3.49 respectively. The speed 
deviation between the both Kaiga generators is plotted in Fig. 3.52 whereas speed deviation 
between Kaiga and is shown in Fig. 3.50. The PSS output at Kaiga-1 is shown in Fig 3.51. 
It can be observed that there is very little oscillation in system parameters. 


Load angle of Machine 1 



Fig 3.15: Load angle oscillation of machine- 1 at Kaiga for reduced NPC system 



Fig.3.16: Power output of machine 1 at Kaiga 
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Fig. 3.22: Speed difference between kaiga-1 and swing machine 
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Fig. 3.27: Speed difference between Kaiga-1 and swing generator 


Load angle of swing generator 



Fig. 3.28: Load angle of swing generator 
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Fig. 3.30: Power output of Kaiga-1 
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Fig. 3.35: Load angle of Kaiga-1, for detailed NPC system with PSS at Kiaga-1 
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Fig. 3.39: PSS output at Kaiga- 
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Fig. 3.40: Terminal voltage at Kaiga- 1 
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Fig. 3.45: Load angle of swing generator 
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Fig. 3.46: Speed difference between Kaiga-1 and infinite generator 
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Fig. 3.47 Load angle of Kaiga-1, for Kaiga system with 16 other generators 



Fig. 3.48: Load angle of swing generator 
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Fig. 3.51: PSS output of Kaiga-1 


PSS output at Kaiga-1 



Fig. 3.52: Speed difference between kaiga-1 and 2 





3.6 Conclusion 


Kaiga generators are working satisfactorily when there is no trip of Kaiga- Kadra and 
Kaiga- Kudasalli lines. But when both of the two lines trip, there is sustained oscillation of 
about 1.2 Hz. This is also verified by simulation studies. When the ratio Ti/T 2 is 0.5 there 
are sustained oscillations. But for the ratio of 0.25, with same loading conditions, the 
oscillations die out rapidly. The PSS tuned is shifting the eigen values to the safe damping 
region there by improving the system dynamic performance. This is also proved through 
non-linear simulation. The determining factors of PSS performance i.e. Ki and K 2 changes 
while the system configuration is altered by incorporating the ICT or adding a new line etc. 
So again an intelligent PSS has to be considered for the system. 
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Chapter 4 


Fuzzy Logic Based PSS Design 

4.1 Introduction 

Conventional PSSs (CPSSs) have been developed using linear control theory to damp 
the oscillation of synchronous machine rotor under transient conditions [10]. Power 
systems are highly non-linear and stochastic in nature because of non-linear loads like 
induction motor, FACTS and HVDC controllers. Also day-by-day increase in demand and 
competition in power market, the power system is becoming very large and complex in 
nature. The control parameters which are optimal for one set of operating conditions may 
not be optimum for another set of operating conditions. This is also demonstrated in 
demonstrated in previous chapter. In some cases, the fixed parameter of PSS may have 
detrimental effect on system dynamic performance. For example when a severe disturbance 
occurs in the system, the AVR forces the field voltage to increase up to the ceiling voltage 
which is known as field forcing [19]. But the CPSS will try to take it out of the ceiling as 
soon as possible, which may lead to system instability. 

Limitations of fixed parameter of conventional PSS has lead to advanced control 
schemes such as self tuning control [22], rule based PSS [14], and fuzzy logic control [14], 
Adaptive control and self tuning techniques require intensive computations and difficult to 
implement. Unlike the classical approach which requires a deep understanding of the 
system, fuzzy logic incorporates an alternative way of thinking, allowing complex system 
to be modeled from accumulated knowledge and experience. A detailed knowledge or 
mathematical model of power system is not required to design Fuzzy Logic PSS (FLPSS). 
The basic feature of fuzzy logic controllers is that the control strategy can be simply 
expressed by a set of rules which describe the behaviour of the controller using linguistic 
terms [22]. Proper control action is then inferred from these rules. In addition, FLC are 
relatively easy to develop and robust as compared to CPSSs. 
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4.2 Fuzzy Logic Controller (FLC) 

In conventional control, the amount of control is determined in relation to a number 
of data inputs using a set of equations to express the entire control process. Expressing 
human experience in the form of a mathematical formula is a very difficult task. Fuzzy 
logic provides a simple tool to interpret this experience into reality. In recent years, fuzzy 
logic has emerged as a powerful tool and is starting to be used in various power system 
applications. The application of fuzzy logic control technique appears to be most suitable 
one whenever a well-defined control objective can not be specified and the system to be 
controlled is a complex one or its exact mathematical model is not available. Fuzzy Logic 
Controllers (FLCs) are robust and have relatively low computational requirements. They 
could be constructed easily using a simple microcomputer. 

Fuzzy logic controllers are rule-based controllers. The structure of the FLC 
resembles a knowledge-based controller except that the FLC utilizes the principle of fuzzy 
set theory [26] in its data representation and its logic. The basic configuration of the FLC 
can be simply represented in four parts as shown in Fig.4. 1 



Fig. 4.1: Schematic diagram of the FLC building blocks 
1. Fuzzification module 

Fuzzification is a process whereby the crisp inputs are mapped onto corresponding 
linguistic (fuzzy) variables. These linguistic variables are used for making decisions. Each 
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linguistic variable or Fuzzy set has a certain membership function. The membership 
function assigns a certain degree of membership or weight to each of the linguistic 
variables. The membership value is nothing but the truth of belongingness of a certain input 
to that fuzzy set or linguistic variable. 

2. Inference mechanism 

Control decisions can now be made based on the linguistic variables. It performs two 
functions. Firstly assigns degree of membership or weight to the variables and secondly 
determines the output decisions or output linguistic variables based on input linguistic 
variables. The definitions of the fuzzy membership functions corresponding to each fuzzy 
set and the necessary rules defining the output goal are included in a knowledge base 

3. Defuzzification mechanism 

This converts the fuzzy set of output variables to crisp numerical outputs. These crisp 
values are the final control inputs to the controlled system. 

4.3 FLC Design 

The design process of an FLC may be divided into the five steps described below. 

(I) Selection of Control Variables 

The selection of control variables (controller inputs and outputs) depends on the 
nature of the controlled system and the desired output. It is more common to use the output 
error (e) and the rate or derivative of the output (e ) as controller inputs. 

(II) Membership Function Definition 

Each of the FLC input signals and output signals, fuzzy variables, have the real line R 
as the universe of discourse. In practice, the universe of discourse is restricted to a 
comparatively small interval [X m j n j, X max j] . The number of fuzzy sets for each fuzzy variable 
varies according to the application. A common and reasonable number of variables are odd 
number. Increasing the number of fuzzy sets results in a corresponding increase in the 
number of rules. A membership function is assigned to each fuzzy set. The membership 
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function maps the crisp values into fuzzy values. A set of membership functions defined for 
seven linguistic variables NB, NM, NS, Z, PS, PM, PB, which stand for Negative Big, 
Negative Medium, Negative Small, Zero, Positive Small, Positive Medium, Positive Big, 
respectively as shown in Fig. 4.2. 



Fig. 4.2: Seven triangular membership functions 

Membership function can be of a variety of shapes, the most usual being triangular, 
trapezoidal, or a bell shaped. The Fig. 4.2 shows a triangular membership function for 
speed deviation as fuzzy variable and is used for design of FLPSS. For simplicity, it is 
assumed that the membership functions are symmetrical and each one overlaps with the 
adjacent functions by 50%. In practice, the membership functions are normalized in the 
interval [-L, L], which is symmetrical around zero. Thus control signal amplitudes (fuzzy 
variables) are expressed in terms of controller parameters (gain). These parameters can be 
defined as: 

r _ 2L 

j ~ X (4.1) 

^ range j 


where X range j defines full range of the control variable Xj, that is, 
rangej ^maxj ~ ^ minj (4.2) 
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Xminj Xmaxj sre the maximum and the minimum values, respectively of the control 
variable Xj. The input and output gains Kj are referred to as FLC parameters. The selection 
of these parameters is usually based on the previous knowledge of the controlled system. 

(Ill) Rule Creation and Inference 

In general, fuzzy systems, as function estimators, map an input fuzzy set to an output 
fuzzy set 8:I n — » I p . Fuzzy rules are the relations between the fuzzy sets. They usually are in 
the form,” if A, then B,” where A is the rule antecedent and B is the rule consequence. Each 
rule defines a fuzzy patch in the Cartesian product A x B (system state space). The 
antecedents of each fuzzy rule describe a fuzzy input region in the state space. This enables 
one to effectively quantize continuous state space so that it covers a finite number of such 
regions [22], In terms of associative memory definition (FAM), each rule represents an 
association (A,; B,). A fuzzy system using two antecedents and one consequence (A, B; C) 
is shown in Fig.3.3. The association (A„ B-,; C,) or the rule of “A,- and B, then C ” maps 
inputs A, B to C\ a partially activated version of C. The corresponding output fuzzy set C 
combines the partially activated sets C 1 , C 2 , C m , that is, 


m 



FAM Rule m 


Fig. 4.3: FAM system architecture for two fuzzy antecedents and one consequence 
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Consider a controller of two fuzzy variables, error (e) and derivative of output (e ), 
and one output fuzzy variable, control strategy (p), each quantized to seven fuzzy sets. This 
leads to a 7 X 7 FAM rule matrix, as shown in Fig. 3.4. Every entry in the matrix represents 
a rule, for example, 

If e is NB and e is NM, then U is NB. 


Error Derivative 










NB 

NB 

NB 
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NB 
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PS 

PM 

PB 

PB 

PB 

PB 


Fig. 4.4: Fuzzy control rule matrix 


The activation of the /th rule consequent is a scalar value w;, which equals the 
minimum of the two antecedent conjucts’ values. For example, if e belongs to NB with a 
membership of 0.3 and e’ belongs to NM with a membership of 0.7, then the rule 
consequence Wj will be 0.3. The knowledge required to generate the fuzzy rules can be 
derived from an off-line simulation, an expert operator, and/or a design engineer. Some of 
the knowledge can be based on understanding of the dynamic system being controlled. 

(IV) Fuzzy Inference 

The well-known inference mechanisms in fuzzy logic are the correlation -minimum- 
encoding and the co-relation-product encoding. The correlation minimum encoding is 
based on the fuzzy outer product notation. The max-min composition operator is denoted 
by the composition operator ‘o\ Thus the fuzzy outer product of the fit row vectors A and 
B, which forms the FAM matrix M: 

M=A t o B, where my =min (ajbj) (4.4) 
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Fig. 4.5: Generation of the output fuzzy set using correlation minimum encoding 

The correlation product encoding uses the standard mathematical outer product of 
the fit vectors A and B to form the FAM matrix M: 

M=A 1 B, where niij= ajbj (4.5) 

(V) Defuzzification Strategy 

Defuzzification is a process of converting the FLC inferred control actions from fuzzy 
values to crisp values. This process depends on the output fuzzy set, which is generated 
from the fired rules [Eq. (4.3)]. The output fuzzy set is formed by either a correlation- 
minimum encoding or the correlation product encoding as discussed in (IV). 

Commonly used defuzzification methods are discussed here: 

( 1 ) Centroid Method 

This procedure also known as center of area or center of gravity method is the most 
prevalent and physically appealing of all the defuzzification methods. It is given by the 
algebraic expression 

* J* jj. ( z).zdz 

J // ( z)dz ^-6) 
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where, jx (z) is the control strategy for the input variable ‘z\ Here J denotes an algebraic 
integration. 

(2) Weighted Average Method 

This method is valid only for symmetrical membership functions. It is given by the 
algebraic expression 

ZjU(z) (4-7) 

where, p (z) is the control strategy for the input variable ‘z\ The weighted average method 
is formed by weighting each membership function in the output by its respective maximum 
membership value. 

(3) Max-membership Principle 

Also known as height method, this scheme is limited to peaked output functions. This 
method is given by the algebraic expression 

jjc 

M(z O) (4.8) 

4.4 FLPSS Design 

A power system stabilizer based on the FLC algorithm can be designed for damping 
the power system oscillations. Since the goal of this application is to stabilize and improve 
the damping of the synchronous machine, speed deviation and active power deviation have 
been selected as controller input [10]. Other options like speed deviation and rate of speed 
deviation, speed deviation and field voltage deviation are also considered. Field voltage 
deviation as one of the inputs is very appealing from the practical point of view, because it 
is very easy to measure as compared to other machine states. As shown in Fig. 4.6 the FLC 
has two input parameters, K w and K p , and one output parameter, Ky. The selection of these 
parameters is usually subjective and requires previous knowledge of the fuzz y control 
variables (input and output signals). Also previous experience of the controlled system 
dynamics is commonly used in the creation of the fuzzy control rules. 
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Fig. 4.6: Schematic diagram of the FLPSS 
4.4.1 FAM table formation 

The previous knowledge of the system parameter variations with conventional PSS is 
used here for the rule base formation. For that, a three phase fault is simulated at the 
generator terminals (the worst case of a fault), and the variation of the required parameters 
is observed. The PSS output is also observed. Same as the FLC design procedures the input 
and output signals of the FLPSS are divided into seven linguistic variables or fuzzy sets, 
such as NB, NM, NS, Z, PS, PM, PB. They have their usual meaning as described above. 
Depending on the variation in conventional case the FAM table was formed, where each 
entry in the table represents a rule. 

1. Rule Base for Aco and AP e as Input Signals (Fig. 4.7) 
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Fig. 4.7: Fuzzy control rule matrix for A© and AP e as PSS inputs 
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2. Rule base for Aco and AE F d as Input Signals (Fig. 4.8) 

Change in Field Voltage 



Fig. 4.8: Fuzzy control rule matrix for Aco and AE F d as PSS inputs 

3. Rule base for Aco and A cb as Input Signals 

When rate change of speed A cb is used as input signal, it can be analyzed as the 
active power signal [2]. So instead of going for observation of the variation in case of the 
fault, the same rule base as rule base 1 is sufficient. Only during implementation, the rate 
change of speed signal has to be multiplied by twice of the inertia constant of the machine 
and inverted before feeding it to the FLPSS. 

4.4.2 FLPSS Parameter Tuning 

Once the proper rules are obtained, the proper parameter ( K w , K p , and KJ) tuning 
should be done in order to achieve good performance. Tuning of FLPSS parameters can be 
a tedious trial-and-error process if not enough information is available about the range of 
the control variables and how they change with different disturbances. 

4.5 Simulink Model for FLPSS 

Fig. 4.9 shows the simulink model used for FLPSS realization. Fuzzy Inference 
System, (FIS) editor of Fuzzy logic toolbox is used for simulation. The inference 
mechanism used is correlation-minimum encoding type (Mamdani type, [22]). So it is 
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called as a Mamdani type FIS system. Centroid method or centre of gravity method is used 
for defuzzification. FLPSS output is given at the AVR summing point in place of CPSS as 
shown in the Fig. 3.10. The inputs to the controller are sampled at a rate of 20ms by a Zero 
Order Hold (ZOH) block. The MUX converts the dual inputs to a vector input to the 
controller. The input and output gains are the PSS parameters. 



Hold 

Fig. 4.9: FLPSS simulink block. 


4.6 Simulation and Results 

The method of defining the fuzzy membership functions used is same as that of FLC 
design, section 4.3. The fuzzy rules in the FAM tables of Fig.4.7 and 4.8 are defined in the 
FIS editor. The input and output ranges are calculated with a three phase fault at generator 
terminal. The different cases as discussed in section 3.5 are simulated with the CPSS 
replaced by FLPSS. 

(a) Reduced NPC system with FLPSS 

The parameters of FLPSS are set as Kp=1.8, K w =1.5, Ku=0.4 obtained by generalized 
GA approach [22]. Load angle of Kaiga-1 and swing generator are shown in Fig. 4.10 and 
Fig. 4.1 1, respectively. Speed deviation between Kaiga-1 and swing generator, and Kaiga-1 
and 2 are shown in Fig. 4.12 and Fig. 4.13 respectively. Power output, terminal voltage 
and PSS output for Kaiga-1 are given in Fig. 4.14-Fig.4.16. From Fig. 4.10 to Fig. 4.16 it is 
observed that, with FLPSS the oscillations are getting damped within less than 5 seconds. 
However, with CPSS, oscillations are damped out in almost 15 second for Tj/T2=0.25. 
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(b) Detailed NPC System 

The FLPSS parameters are maintained at the same value, and the simulation was 
performed for the detailed model. Load angles of Kaiga, swing generator and Kadra 
generator are shown in Fig. 4.17, Fig. 4.18 and Fig. 4.20, respectively. Speed difference 
between Kaiga-1 and swing generator and PSS output at Kaiga-1 are shown in Fig. 4.19 
and Fig. 4.21 respectively. FLPSS also improves the system performance in this case as 
compared to CPSS. 

c) Kaiga System with Other 16 Generators 

Other generators remain the same; FLPSS is installed only at Kaiga-1. The load 
angles of Kaiga-1, swing generator and Kadra hydro are shown in Fig. 4.22, Fig. 4.23 and 
Fig 4.24, respectively. Speed deviation between Kaiga-1 and Kadra, Kaiga-1 and 2 
generators are shown in Fig. 4.25 and Fig. 4.26 respectively. Kiaga-1 terminal voltage is 
given in Fig. 4.27. The Fig. 4.22- Fig. 4.27 shows that system FLPSS is able to damp the 
system oscillations very quickly following a three phase fault and subsequent opening of 
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Fig. 4.10: Load angle of Kaiga-1 for reduced NPC system 
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Fig. 4.12: Speed difference between Kaiga-1 and swing generator 
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Fig. 4.14: Power output of Kaiga-1 
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Fig. 4.20: Load angle at Kadra generator 
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Fig. 4.22: Load angle at Kaiga-1 with fuzzy PSS for other 16 gen connected 


92 





Angle in degrees 






leed in PU— 


speed deviation between Kaiga and Kadra generator 


0.09 
0.08 
0.07 
a 0-06 
0.05 
0.04 
S' 0.03 
0.02 
0.01 
0 



25 30 35 40 45 50 

Time in seconds — > 


55 60 


Fig. 4.25: Speed deviation between Kaiga-1 and Kadra generator 
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Fig. 4.26: Speed deviation between kaiga-1 and 2 






Fig. 4.27 terminal voltage of Kaiga-1 


4.7 Conclusion 

FLPSS can be used as a better alternative for CPSS to damp the rotor oscillations. 
The FLPSS input and output parameters do not vary for a change in system configuration. 
The range of the control inputs change with change in system configuration, but to certain 
extend the same PSS parameters are able to damp the rotor oscillations. Only for extreme 
contingency cases, the PSS may not be able to damp the oscillations within reasonable time 
if the input parameters range exceeds the specified range. For this reason adaptive FLPSS is 
preferred, where the input and output parameters are determined depending on the severity 
of contingency. The input and output parameters are optimized for a three phase fault at 
generator terminal of reduced NPC system. 
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Conclusions 


Chapter 5 


5.1 General 

Because of long transmission lines, operation of the generators near to their full load 
capacity and poor tuning of AVR parameters, present day practical power systems are 
facing the problem due to poor damping of small signal oscillations. Load centers and 
generating stations are placed far from each other because of various reasons. So long 
transmission line (weak system) is inevitable. Because of heavy demand mostly in 
developing countries, the generator has to operated near the full load capacity. Higher 
settings of AVR increases the field forcing capability, thereby increases the steady state 
stability limit. But simultaneously magnifies the oscillations by reducing the damping 
torque [15]-[16], The problem of long transmission line can be solved by application of 
FACTS devices; which is beyond the scope of the present thesis. This thesis identifies poor 
tuning of AVR parameters as the main cause of undamped oscillations in a practical system 
and prescribes tuning of PSS for damping these oscillations. 

In this thesis, the causes of small signal oscillations in Kaiga generators were found 
out. The remedial measure to damp these oscillations was provided by tuning the PSS and a 
Fuzzy Logic PSS is proposed to replace the conventional PSS. 

The main findings of the studies carried out in chapter-2 are as follows: 

• In the 9 bus system, all the eigen values had negative real parts but the damping 
was poor with one critical mode in the system. Application of strip eigen value 
placement method and pseudo decentralization method improved the system 
damping. 

• In the 39-bus system there were three under damped modes with damping ratio 
< 0.05. Application of pseudo decentralization method improved the system 
damping ratio. But strip eigen vale method could not be applied to this system. 
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In chapter 3, the SMIB approach was used for tuning of Kaiga PSS, as there is no lead-lag 
block as there in practice. The main findings of this chapter are 

• The PSS signal had to overcome considerable phase lag in the GEP(s) block to 
produce the final torque on the rotor in phase with speed deviation signal. 

• Eigen value analysis showed that there are two critical modes and both the 
generators had maximum participation in AS and Aco states corresponding to these 
modes. After introduction of PSS the damping ratio had improved to safe values 

• With Tj/T 2 setting of 0.5, the simulation of the reduced NPC system exhibited 
similar oscillations at 1.2 Hz. But with Tj/T2=0.5, the oscillations damped out. This 
proved the improper tuning of AVR at Kaiga. After inclusion of PSS with the 
settings found through eigen value analysis, the oscillations were damped out. 

• The detailed model considering the dynamics of near by hydro generators had 
similar oscillating response without PSS. After inclusion of PSS the oscillations 
were damped out. 

• The NPC system was tested for a three phase fault case, with and without PSS. The 
PSS was able to damp the oscillations. Which proved that, this setting of PSS can 
also improve the system performance during severe fault conditions? 

Chapter-4, CPSS was replaced by FLPSS. The main findings of this chapter are 

• The rule base was prepared for the FLPSS through the simulation of the NPC 
system with a 3 phase fault at the generator terminal. So a detailed knowledge of the 
system dynamics is not required for finding the rule base, as it can also be designed 
from operating experience of the system 

• The FLPSS was able to damp the oscillations within few seconds, so it can be used 
as an alternative to CPSS. 
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5.2 Future Scope of Research 


During the course of research, following areas were found out for further research and 
development 

> The decentralization method can be improved further, so that an effective method to 
find out the Q matrix can be obtained. This will reduce the burden of arbitrary 
selection and variation of PSS parameters. So that practical application of this 
method will be feasible. 

> The actual settings of nearby generators to Kaiga were not available for analysis. 
The analysis can be made more realistic if those settings are obtained from the site. 

> FLPSS needs some adaptive technique, so that it can take care of different operating 
conditions. This remains as a future area of research. 
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Appendix - A 


WSSC 3 -Machine, 9-Bus System 
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Fig A-l one line diagram of WSCC 9-bus system. 
Table A-l Base case load flow Results 


Bus 

Voltage 

Phase 

(deg) 

Pgen 

Qgen 

Pload 

Qload 

1 

1.040 

0.00 

0.7164 

0.2705 

0.0000 

0.0000 

2 

1.025 

9.28 

1.6300 

0.0665 

0.0000 

0.0000 

3 

1.025 

4.66 

0.8500 

-0.1086 


Hiltltltfsii 

4 


-2.22 

0.0000 

0.0000 


MW 

5 


-3.99 

0.0000 

0.0000 



6 

WBESSKM 

3.72 

0.0000 


■sEIiIsIiM 


7 

■BS9| 

3.72 


i 



8 

1.016 

0.73 


1 I 



9 

1.032 

1.97 


| 
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Table A-2 Line Data 


From 

To 

R_ line 

X line 

Y charging 

Tapping 

4 

1 

0.0000 

0.0576 

0.0000 

1 

7 

2 

0.0000 

0.0625 


1 

9 

3 


■OHS 

0.0000 

1 

4 

6 



0.1580 

0 

4 

5 




0 

5 

7 




0 

6 

9 


0.1700 


0 

7 

8 




0 

8 

9 


jpfflKPjP 


0 


Table A-3 Machine Data 


Gen 

H 

x d 

X d ’ 

x q 

X q ’ 

T d0 ’ 

V 

l 

23.6 

0.1460 

0.0608 

0.0969 

0.0969 

8.96 

0.31 

2 

6.4 

0.8958 

0.1198 

0.8645 

0.1969 

6.00 

0.54 

3 

3.0 

1.3125 

0.1813 

1.2578 

0.2500 

5.89 

0.60 


Table A-4 Exciter Data 

The IEEE type-1 exciter having identical parameters has been used in all the three 
generators. The exciter data is given below: 

K a = 20.00 
T a =0.2 sec 
K g = 1.0 
T e =0.314 sec 
K f =0.063 
T f =0.35 sec 
S E (E fd ) =0.003 9e 1 555E F< i 
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0.000 
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Table B-2 Line Data 


From 

To 

R_line 

X_line 

Y_charging 

Tapping 

39 

30 

0.0007 

0.0138 

0.0000 

1 

39 

5 

0.0007 

0.0142 

0.0000 

1 

32 

33 

0.0016 

0.0435 

0.0000 

1 

32 

31 

0.0016 

0.0435 

0.0000 

1 

30 
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0.0180 

0.0000 
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29 
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0.0156 

0.0000 
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25 
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1 

23 
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0.0000 

1 

22 
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0.0143 

0.0000 

1 

20 
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0.0000 
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0.0000 

1 

16 
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0.0000 

1 

37 
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0.0013 

0.0173 

0.3216 

0 

37 
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0.1319 

0 

36 
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0 

36 

21 
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0.2548 
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Tabic B-3 Machine data 


Gen 

I-I 

x d 

x d ’ 

x« 

Xa* 

Tdo’ 

T ’ 

1 qo 

1 

30.3 

0.2950 

0.0647 


0.0647 

6.56 

1.50 

2 

500.0 

0.0200 

0.0060 


0.0060 

6.00 

0.70 

3 

35.8 

0.2495 

0.0531 


0.0531 

5.77 

1.5 

4 


0.3300 

0.0660 


0.0660 

5.40 

■EEI 
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0.2620 

0.0436 
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0.0436 
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0.0500 


0.0500 
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0.0490 


0.0490 

5.66 
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24.3 

0.2900 

0.0570 


0.0570 

6.70 

0.41 
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liWfifl 

0.0570 


0.0570 

4.79 
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0.0040 


0.0040 
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For all machines R a = 0 


Table B-4 AYR data 


Machine 

no. 

Ka 

T a 

K e 

T e 

K f 

Tf 

1 

5.0 

0.06 

- 0.05 

0.25 

0.04 

1.000 

2 

6.2 

0.05 

0.63 

0.41 

0.06 

0.500 

3 

5.6 

0.06 

- 0.02 

0.5 

0.08 

1.000 

4 

5.0 

0.06 

- 0.05 

0.5 

0.08 

1.000 

5 

40.0 

0.02 

- 0.04 

0.785 

0.03 

1.000 

6 

5.0 

0.02 

1.00 

0.471 

0.08 

1.250 

7 

40.0 

0.02 

1.00 

0.730 

0.03 

1.000 

8 

5.0 

0.02 

- 0.05 

0.528 

0.09 

1.260 

9 

40.0 

0.02 

1.0 

1.4 

0.03 

1.000 

10 

5.0 

0.02 

1.0 

1.4 

0.03 

1.000 
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Appendix - C 


System and Machine Parameters for Kaiga NPP 


1 . M/C inertia constant (H) 191 MJ/MVA 

2. M/C Voltage 16 5 kv 

Table C-l Machine reactances in PU on machine voltage and power base 



Unsaturated 

Saturated 

Direct axis synchronous reactance (X d ) 

1.94 

1.72 

Quadrature axis synch. Reactance (X q ) 

1.84 

1.63 

Direct axis transient reactance (X d) 

0.28 

0.238 

Quadrature axis transient reactance (X q ) 

0.82 

0.74 

Direct axis subtransient reactance(X d) 

0.209 

0.167 

Quadrature axis subtransient reactance(X q ) 

0.23 

0.184 

Potier reactance (Xp) 

0.209 

0.167 


Table C-2 Machine time constants in second 
Open circuit Direct axis transient time constant (T do) 6.3 

Open circuit Quadrature axis transient time constant (T qo ) 2.0 

Open circuit direct axis subtransient time constant (T do) 0-04 

Open circuit Quadrature axis subtransient time constant (T qo ) 0.17 


Table C-3 Exciter and AYR parameters 


V 0 =steady state gain 
Vp=Proportional gain 

y,„ f 

T R := transducer time constant 
T a =AVR time constant 
Tb=AVR time constant 


Range 

Typical values 

10-126 

40 

0.8-25.6 

10 

8-257 

200 

- 

0.025sec 

- 

1 sec 

- 

10msec 
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Ts=Thyristor time constant 


3.3(j.sec 


Table C-4 Line and Transformer data 

Line data in ohms/km/circuit, in 100MVA and 220kV base 



R, 

X, 

Ro 

X 0 

For single circuit lines 

0.08 

0.399 

0 

0 

For all double circuit lines 

0.029792 

0.322 

0.16192 

1.24 

Transformer data on Machine base 

X, 

X 0 

Xi/R, 

X 0 /Ro 

All two winding transformers 

0.14 

0.14 

20 

20 

Three winding transformer on machine base 

X, 

x 2 

X 0 

Primary-secondary reactance 

45% 

45% 

95% 

Secondary to tertiary reactance 

12.5% 

12.5% 

95% 

Primary to tertiary reactance 

30% 

30% 

95% 


110 



